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ABSTRACT 

This  report  comprises  two  memoranda,  each  with  its  own  abstract.  The 
reader  is  referred  to  these. 


PREFATORY  NOTE 


In  the  course  of  research  under  this  contract  on  methods  for  more  efficiently 
utilizing  antennas  in  fading  radio  communications,  such  as  typified  by  tropospheric 
scatter,  it  was  desired  for  comparison  purposes  to  carry  out  detailed  performance 
estimates  for  angle  diversity  systems.  Such  systems  have  been  proposed  on  the 
basis  of  experimental  observations  that  the  apparent  angle  of  arrival  of  such  fading 
signals  tends  to  wander  over  several  beamwidths  at  the  receiving  point;  it  then 
appears  economical  in  equipment  to  consider  a  diversity  system  based  on  the  use 
of  a  single  large  reflector  with  a  number  of  feeds  arranged  to  produce  squinted  beams . 
Each  feed  is  connected  in  such  a  system  to  a  separate  receiver,  and  the  outputs  of 
the  receivers  are  combined  in  the  appropriate  diversity  combining  circuitry,  for 
exarrple,  so-called  maximal- ratio  combining. 

A  number  of  analyses  of  angle  diversity  have  been  reported  to  date.^^”^^ 


S.  Stein,  D.  E.  Johansen,  andA.  W.  Starr,  ''Theory  of  Antenna  Performance  in 
Scatter-Type  Reception".  Hermes  Electronics  Co.  Report  M783,  AFCRC-TR-59- 
191.  30  Sept.  1959  (Appendix  F). 

2.  A.  B.  Crawford,  D.  C.  Hogg,  and  W.H.  Kummer,  "Studies  in  Tropospheric 
Propagation  Beyond  the  Horizon",  BSTJ,  38,  pp.  1067-1178  (1959). 

3.  J.  H.  Vogelman,  J.  L.  Ryerson,  and  M.  H.  Bickelhaupt,  "Tropospheric  Scatter 
System  Using  Angle  Diversity",  Proc.  IRE,  47,  pp.  688-696  (1959). 

4.  R,  Bolgiano,  Jr.,  N.H.  Bryant,  and  W.  E.  Gordon,  "Diversity  Reception  in 
Scatter  Communications  with  Emphasis  on  Angle  Diversity",  Cornell  Univ. , 

Final  Report,  Part  I,  on  AF30(602) - 17  17,  Jan.  1958,  . 

5.  C,  Chu,  A,  H,  Wren,  and  J.  LaRue,  "Evaluation  of  the  Pincushion  System", 

Univ.  Mich,  Rad,  Lab.  Rept.  2872-1-T,  RADC-TN-6O- 50,  Feb.  I960, 

6.  W.R.  Richard  and  M.  H.  Bickelhaupt,  "Multiple  Angle  Diversity  Design 
Considerations",  Rome  Air  Development  Center,  RADC-TN-60-22,  April  I96O,. 
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Problems  which  are  discussed  in  this  scientific  report  were  first  noted  in  an 
earlier  analysis^ where  it  was  disturbing  to  find  that  casual  angle  diversity 
calculations  led  to  results  which  were  obviously  inconsistent  with  conservation 
of  energy.  The  same  type  of  mathematics  appears,  or  seems  to  be  implied,  in  all 
the  other  discussions  of  angle  diversity.  The  problem  which  arises  is  the 

following:  The  "casual  diversity  calculation",  referred  to  above,  is  one  which 
assumes  that  each  of  the  various  beams  of  the  receiving  array  collects  energy 
from  the  incident  field,  essentially  independently  of  the  presence  of  the  other 
beams.  That  is,  one  calculates  the  received  power  in  each  beam  from  a  knowledge 
of  the  incident  field,  or  equivalently  of  the  sources  creating  the  field,  by  applying 
the  reciprocity  theorem  with  the  radiation  pattern  of  each  beam  only,  assuming 
that  there  are  no  interaction  effects  due  to  the  presence  of  the  other  beams  which 
would  diminish  the  estimate  of  received  power.  The  basis  for  this  assumption  is 
that  the  beams  do  not  greatly  overlap.  However,  althoagh  for  reasons  of  physical 
construction,  standard  multiple  feed  microwave  antenna  systems  usually  cannot 
produce  beams  which  cross  over  at  levels  higher  than  about  their  3  db  points,  it  is 
not  at  all  clear  that  with  some  ingenuity  useful  beams  could  not  be  designed  which 
would  overlap  even  more  greatly.  Surprisingly,  the  analysis  in  this  report  indicates 
that  even  for  beams  crossing  over  at  their  3  db  points  the  modifications  indicated  by 
the  analysis  can  be  quite  significant. 

Of  course,  somewhat  paradoxically,  all  the  angle  diversity  analyses  have 
shown  an  awareness  of  the  need  for  taking  account  of  beam  overlap  insofar  as 
estimating  the  degree  of  correlation  of  the  fluctuations  in  the  energies  received 
on  the  various  beams,  in  the  case  of  fading  propagation  channels  typified  by 
tropospheric  transhorizon  ("scatter")  systems.  However,  the  issue  to  be  discussed 
in  this  report  is,  in  fact,  quite  separate  from  the  fluctuating  or  nonfluctuating 
nature  of  the  incident  field.  It  is  rather  a  statement  of  inherent  limitation  in  the 
character  of  the  multiple  beam  antenna  systems,  due  to  the  necessity  of  the 
conservation  of  energy  principle.  In  short,  the  relations  to  be  described  here 
refer  to  a  limitation  on  the  energy  which  can  be  collected  at  each  instant  by  the 
antenna  system  in  a  fading  environment,  and  do  not  refer  at  all  to  time  fluctuations 
of  such  collected  energy. 

The  nature  of  the  problem  under  discussion  becomes  obvious  by  considering 
as  an  example  the  extreme  limiting  case  of  completely  overlapping  beams.  For 
example,  suppose  one  has  an  antenna  system  with  one  feed,  but  chooses  to  put  M 
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different  taps  on  the  feedline  at  various  points  along  the  circuit.  Looking  at 
everything  between  those  taps  and  the  radiation  patterns  as  a  "black  box  antenna 
system",  which  one  is  certainly  free  to  do  (if  only  to  be  different),  gives  an  antenna 
system  with  M  "feed-lines"  entering  the  "black  box",  each  exciting  a  perfectly  good 
radiation  pattern  when  excited,  with  all  the  other  taps  terminated  in  matched  loads 
(which  is,  indeed,  how  one  defines  the  beam  patterns  in  any  multiple  beam  system). 
Thus,  we  would  observe  that  our  black  box  is  a  multiple-beam  antenna  system  which 
has  M  identical,  completely  overlapping  beams.  It  is  certainly  clear  on  the  one  hand 
that  one  cannot  buy^  something  for  nothing,  that  in  reality  there  exists  only  the  one 
radiation  pattern;  that  if  in  receiving,  the  outputs  of  the  taps  are  combined  in  an 
optimum  linear  combining  network,  the  net  results  will  at  best  be  the  same  as  if  one 
had  initially  adopted  the  more  conventional  attitude  that  there  exists  only  the  one 
feed,  and  therefore,  only  one  real  beam  collecting  energy  and  one  terminal  point. 

In  fact,  of  course,  if  one  excites  one  of  the  taps  in  our  example,  some  fraction  of 
the  energy  will  disappear  into  the  loads  associated  with  the  other  taps  and  only  the 
remaining  fraction  into  the  radiated  beam.  In  other  words,  the  presence  of  these 
other  taps,  or  alternatively  the  fact  that  the  beams  overlap  as  much  as  they  do, 
creates  an  insertion  loss  between  each  tap  and  its  associated  beam.  It  is  this 
insertion  loss,  which  appears  in  reception  as  well  as  in  transmission,  which  is 
inherent  in  the  principle  of  conservation  of  energy,  and  for  which  we  seek  a 
characterization. 

In  the  extreme  example  just  cited,  the  problem  was  obvious  and  the  correction 
factors  equally  obvious.  The  challenge  which  appeared,  however,  since  we  wanted 
to  analyze  angle  diversity  configurations  in  which  the  beams  crossed  over  at  varying 
squint  angles,  and  hence  with  varying  degrees  of  overlap,  was  to  produce  a 
mathematical  formulation  and  description  of  such  a  system  which  will  automatically 
and  routinely  account  correctly  for  these  limitations.  Such  a  formulation  would  then 
be  completely  dependable  for  performance  predictions  for  angle  diversity  systems, 
regardless  of  how  peculiar  the  feeding  system  (or  the  way  in  which  the  inventor 
chooses  to  regard  the  system). 

However,  this  mathematical  formulation  appeared  possibly  to  be  of  much  more 
general  significance  to  designers  of  multiple  beam  antenna  systems.  For  this  reason, 
while  the  analysis  performed  under  this  contract  achieved  its  goal  once  a  formulation 
giving  the  correct  results  for  angle  diversity  had  been  achieved,  it  seemed  desirable 
to  further  explore  these  results  in  the  more  general  context.  Although  this  further 
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exploration  was  carried  out  external  to  the  contract,  for  the  sake  of  completeness 
the  entire  contents  are  being  reported  in  this  scientific  report,  with  reference  to 
the  more  general  context. 

This  scientific  report  consists  of  two  engineering  memoranda  which  are 
entitled: 

Part  I;  Radiation  Efficiencies  (with  Application  to  Angle  Diversity) 

Part  II:  Characteristics  of  the  Cross -Coupling  Matrices 

The  contents  of  each  are  perhaps  obvious  from  the  titles.  It  is  our  belief,  as 
stated  already,  that  a  significant  new  insight  has  been  obtained  through  this  work 
into  many  of  the  problems  of  designing  multiple  beam  antenna  systems.  Never¬ 
theless,  as  will  be  apparent  after  reading  this  report,  much  further  useful  work 
remains  to  be  done  if  the  full  significance  of  these  results  is  to  be  obtained  with 
respect  to  particular  system  applications.  These  conclusions  are  reiterated 
perhaps  more  meaningfully  at  the  ends  of  the  two  memoranda. 
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PART  I  RADIATION  EFFICIENCIES 
(With  Applications  to  Angle  Diversity) 
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ABSTRACT 

A  general  multiple-beeun  antenna  syatem  is  considered,  in  vhich  the  beam 
patterns  overlap.  The  degree  of  overlap  is  defined  by  a  cross -correlation  type 
of  integral  which  inclvides  phase  factors  as  well  as  polar  diagrams  in  the  des¬ 
cription  of  the  beams.  It  is  shown  that  when  beam  overlap  exists,  conservation 
of  energy  in^jlies  the  existence  of  unavoidable  cross-couplings  between  the  feed¬ 
lines,  and  a  related  limitation  on  the  reidlatlon  efficiency  (or  corresponding 
receiving  cross  section)  of  any  single  beam  when  excited  alone.  Certain  problems 
relating  to  the  angle  diversity  concept  are  discussed,  by  way  of  exemple.  In 
addition,  some  unsolved  corollary  problems  with  possibly  important  design  applica¬ 
tions  are  outlined. 

In  Part  II  (a  ccmpanion  paper),  the  scattering  matrix  comprised  of  the  array 
of  cross-coupling  factors  between  the  various  feed-lines  is  considered.  Some  pro¬ 
perties  are  derived,  and  Isplicatlons  discussed.  In  addition,  further  unsolved 
problems  relating  to  characterization  of  the  scattering  matrix  are  also  outlined. 


1 .  Introduction : 
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There  are  many  multiple-beeun  antenna  systems  in  use  or  contemplated  in  which  the 
radiating  structure  may  be  regarded  as  driven  by  a  number  of  feed- lines,  each  feed¬ 
line  corresponding  to  excitation  of  a  particular  beam.  The  various  beams  may  be  exci¬ 
ted  simultaneoixsly  or  sequentially,  the  assunqption  being  made  below  that  the  radiating 
structure  is  not  changed  physically  sinqjly  because  a  particular  feed-line  is  not  exci¬ 
ted.  When  there  is  overlap  among  any  of  the  beams  which  the  structure  is  designed  to 
radiate  (overlap  in  a  sense  to  be  defined  below),  it  becomes  apparent  that  the  corres¬ 
ponding  excited  portions  of  the  radiating  structure  are  likely  to  overlap,  and  there¬ 
fore  that  crcss-coiq)ling  may  exist  between  the  corresponding  feed- lines.  In  fact,  it 
will  be  shown  below  that  such  beam  overlap  renders  cross-csoupling  unavoidable,  and  it 
is  the  object  of  the  paper  to  describe  various  iD5>lications  in  terms  of  the  beam  pat¬ 
terns.  Most  of  the  attention  in  this  Part  I  is  given  to  deriving  the  basic  equations, 
and  subsequently  discxxssing  radiation  efficiency  problems,  with  particular  attention 
to  the  angle  diversity  concept.  In  Part  II^^^,  the  detailed  characterization  of  the 
cross -couplings  is  treated. 

2.  Scattering  Matrix  Representation  of  the  Antenna  System. 

Consider  an  antenna  system  con5)rising  M  feed-lines,  feeding  into  a  radiating 
structure  which  we  regard  as  a  junction.  We  define  (Fig.  l)  the  k  beam  pattern 
as  that  beam  radiated  when  a  generator  is  exciting  the  k  feed-line,  and  it  and 
nil  other  feed-lines  are  terminated  in  matched  loads  which  absorb  any  energy  reflected 
back  from  the  Junction. 

Rather  than  describe  the  antenna  system  in  terms  of  voltages  and  currents,  or 
iii5>edances ,  it  will  be  useful  to  eitploy  an  incident  and  reflected  wave  representa¬ 
tion  (similar  to  the  description  of  transmission  lines  in  terms  of  forward-  and 
backward-traveling  waves).  Specifically,  let  us  consider  a  steady-state,  single 

Stein  and  J.E.  Storer,  "On  Cross -Coupling  In  Multiple  Beam  Anteniwis,  Paurt  II. 

Characteristics  of  the  Cross-Coiqpling  Matrices,"  AEM-238,  2  March  I96I. 


-1- 


ARM-237 

frequency  situation  (harmonic  time  dependence).  We  will  denote  the  an5>lltudes  of 
the  waves  traveling  towards  the  Junction  in  each  of  the  feed-lines  by  the  set  of 
conqjlex- valued  numbers  M  and  the  waves  traveling  away  from  the 

Junction  at  the  same  time  by  the  set  ,  k  =  1,  M.  We  define  the  scales 

of  the  Individual  and  y^^  in  each  line  such  that  a  unit  magnitude  corresponds  to 
a  unit  power  (l  watt)  level  in  the  wave  (power  transfer  by  the  wave  across  a  cross- 
section  of  the  line).  The  phase  factors  of  the  x^^  and  y^^  refer  to  scxne  arbltrsury 
reference  point  in  each  line. 

We  further  ass\me  the  antenna  system  contains  no  non-linear  elements.  Each 
of  the  y^  is  therefore  related  to  the  Xj^  by  some  linear  sum,  of  the  form 


M 

^  ®km  ^m 
m=l 


(1) 


where  the  are  so-called  scattering  coefficients.  If  we  write  the  arrays  of  the 


Xj^  and  yj^  as  single -column  matrices  (vectors), 


X  »- 


(2) 


and  denote  the  M  x  M  array  of  the  Sj^  as  the  scattering  matrix,  S  (S^^  being  the 
element  in  the  k^^  row  and  column),  we  can  rewrite  Eq.  1  as  the  matrix  eqxiatlon 


Sx 


(3) 


Next,  let  us  define  the  M  beams  in  detail  as  follows:  Consider  a  generator  In 
the  k  feed-line,  exciting  a  unit  anqplltude  wave  incident  on  the  radiating  structure 


-2- 


ARM-23T 


in  the  feed-line,  and  with  all  other  feed-lines  terminated  in  matched  loads 
so  that  there  is  no  mechanism  to  produce  incident  waves  in  these  other  feed-lines, 
i.e., 

Xj  =  0  ,  J  k 

Under  these  circimsteuices,  let  the  conqplex- valued  field  radiated  by  the  antenna 
structure  be  described  (in  the  far-zone)  by 


exp(i  ^  r) 


(4) 


The  vector  Q  represents  the  spatial  angle  coordinates,  and  r  the  radial  distance 

from  a  reference  origin  located  at  the  antenna  structure.  The  explicit  e:iq)(i  r) 

dependence  on  r  is  the  usiial  outgoing  spherical  wave  form  describing  the  f„r-zone 

field,  with  the  form  of  the  angular  dependence  described  by  The  vector  form 

for  and  denotes  the  possibility  of  two  crossed  linear  polarization  components 

(e.g.,  "horizontal"  and  "vertical")  at  each  point  in  the  far-zone.  (Each  polariza- 

tion  con^jonent  is  also  con5)lex-valued. )  We  label  the  angular  dependence  ^(®)  El’S 
th 

the  k  beam  pattern .  We  will  assume  R,  (o)  to  be  conveniently  normalized  such  that 


'^(0)  •  Rj^(0)  dO  =  1 


(5) 


where  is  the  "admittance  of  free-space";  the  asterisk  denotes  a  conqplex 

conjixgate;  the  usual  vector  dot  product  (scalar  product)  is  indicated;  and  the  inte¬ 
gration,  as  indicated  symbolically,  is  over  all  spatial  angles.  With  this  normali¬ 
zation,  the  magnitude  of  the  coefficient,  q  ,  is  the  relative  an5>litude  of  the  field, 

iC 
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and  its  phase  refers  the  phase  of  the  field  to  that  of  the  vave  incident  in  the 
feed-line  and  to  the  definition  of  a  phase  reference  in  Rj^(o).  In  particular,  the 
total  radiated  power  in  the  far-zone  field  of  Eq.  (4)  is  given  by 


/O  -Mi'  ^  -s  uk  ^ 

r^  Ej^(O)  •  Ej^(o)  d© 


(6) 


Using  Eq.  (4),  and  the  normalization  in  Eq.  (5),  this  reduces  to 
w  *  I  '2 


(7) 


where  bars  indicate  the  magnitude  of  complex  quantities.  Since  unit  power  is  ihci- 

I  |2 

dent  upon  the  Jiinction,  qj^  may  be  considered  to  represent  the  radiation  efficiency 
th 

for  this  k  beam.  Also,  with  the  particular  excitation  euad  loading  specified, 

12 
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most  represent  energy  carried  away  from  the  junction  in  reflected  waves 


in  the  several  feed-lines  and  absorbed  in  their  terminating  matched  loads;  alter¬ 
natively  described,  insertion  loss  between  the  input  wave  on  the  k^^ 

th 

feed- line,  and  the  radiated  k  beam. 


3.  Beam-Coupling  Factors 

Having  defined  the  M  beam  patterns,  we  now  define  a  measure  of  the  beam  over¬ 
laps.  The  appropriate  measure,  as  will  be  seen  later,  is  the  set  of  beam-coupling 
factors  defined  by 

^kj  -  T  "  (8) 

Clearly, 

^kj  “  ^Jk 
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and  by  Eq.  (5), 


6  =1 


kj 


<  1 


(9b) 

(9c) 


It  is  ingjortant  to  note  that  Eq.  (8)  does  not  define  beam  enupTlwg  simply  in  terms 
of  overlap  of  the  usxial  polar  diagrams.  The  fact  that  the  Rj^(6)  are  complex- valued 
can  be  extreuely  important,  since  rapid  variations  in  phase  in  the  product  being 
integrated  can  cause  the  integral  to  vanish.  Such  rapid  variations  can  occur,  for 
exanple,  if  the  locations  on  the  radiating  structure  of  the  major  cvirrents  exciting 
the  two  beams  are  sufficiently  disjoint.  For  instance,  if  a  set  of  currents  gives 
rise  to  a  beam  R(0),  and  this  set  of  currents  is  translated  by  a  distance  d,  the  new 
pattern  (?  and  r  retaining  their  original  meaning)  will  have  the  form 


esq) 


[■ 


i2jt  ^  cos 


where  ♦  is  the  angle  between  the  direction  of  the  observing  point  in  space,  and  the 
direction  in  which  the  currents  were  translated.  As  another  exaaple,  if  we  (l6osely) 
use  the  Fourier  Transform  relationship  between  the  beam  pattern  and  the  aperture 
plane  fields  for  a  microwave  "aperture-type"  anteima,  the  integral  In  Eq.  (8)  is 
readily  shown  to  be  equivalent  to  an  Integral 


/ 


F,(s) 


Fj(s)  ds 


where  s  represents  coordinates  in  the  apert\ire  plane,  and  F(s)  the  aforementioned 
apert\ire  fields.  But  then,  if  Fj^(8)  and  Fj(s)  are  respectively  non-vanishing  only 
over  disjoint  sections  of  the  aperture  plane,  the  Integrsd  vanishes.  In  actual 
fact,  such  disjointness  may  be  more  a  design  ideal  than  a  realized  flact,  due  to 
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induced  currents,  "leakages",  etc;  and  many  designs,  e.g.,  phased  arrays,  of 
course  do  not  have  such  disjointness  to  begin  with,  (in  fact,  beam  pattern 
measiirements  may  represent  a  more  effective  engineering  approach  to  determining 
the  actual  extent  of  beam-couplings  than  trying  to  determine  what  actually  are  the 
totality  of  cvirrents  associated  with  any  given  beam  pattern.) 

Any  pair  of  beams  for  which  the  corresponding  beam  coupling  factor  vanishes 
will  be  said  to  be  decoupled. 


4.  Conservation  of  Energy  Equations 


We  can  now  proceed  to  the  main  results  of  this  paper.  Consider  again  the 
steady-state  situation  envisaged  in  Section  2,  with  now  a  general,  excitation  of 
the  feed-lines  such  that  waves  are  incident  on  the  radiating  structure  in  all  feed¬ 
lines,  with  the  set  of  an^ditudes  given  by  the  column  matrix,  x.  As  described  by 
Eq,  (3),  the  associated  reflected  waves  in  all  feeding  lines  are  given  by  the 
column  matrix,  y  =  Sx.  Cleeirly,  however,  we  can  now  complete  the  picture;  Corres- 
ponding  to  the  incident  wave  of  conqplex  amplitude  in  the  k  line,  there  is  a 

radiated  field, 

-k  -»  -*  -»  exp(i  ^  r) 

^  ‘Ik 


and  hence,  by  superposition,  the  total  radiated  field  associated  with  x  is 


e3q>(i  Y 
r 


(10) 


Further,  the  power  flowing  towards  the  radiating  structure  is  the  svm  of  the  powers 
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(11) 


In  matrix  notation,  we  define  the  transpose  of  the  column  matrix  x  hy  the  row 
T 

matrix  x  , 


T 

X 


(Xi,  Xg 


(12) 


and  by  the  usual  laws  of  matrix  multiplication. 


TJ  T* 

P .  ,  ,  =  X  X 

total 


(13) 


Similarly,  the  total  power  carried  by  the  reflected  waves  in  the  feed-lines  is 
simply, 

T. 


(111) 


Now,  we  can  also  define  the  transpose  of  S, 


(15) 


And  finally,  let  us  now  define  the  conjugate  transpose  of  S,  denoted  by  the  matrix. 


+ 

S  ,  as 


s  =  (s^)*  -  (s  r 


(i&) 


and  similarly  define  the  conjugate  transpose  column  vectors 


+  f  Tx* 

X  =  (x  ) 


+  ,  Tx< 

y  =  (y  ) 


(l6b) 
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We  also  recall  that  for  two  matrices,  A  said  B, 


Thus  we  may  rewrite  Eqs.  (13)  and  (lU)  as 


P.  .  ,  «  X  X 
total 


’'refl 


and  recalling  Eq.  (3), 


^refl  ■ 


+  _+_ 

X  S  Sx 


Finally,  the  total  radiated  power,  using  Eqs.  (lO),  (9)  and  (8),  is 


r^  ^  E  (0)  •  E(0)  dO 


M 


^  »  » 

k,J=l  J  ■ 


i 


(0)  .  Rj(0)  dO 


If  we  further  define  an  M  x  M  matrix,  P,  whose  elements  sure 


we  can  write 


P  =  x^  r  X 
rad 


(17) 

(18) 

(19) 


(20) 


(21) 


(22) 


(23) 
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The  lav  of  conservation  of  energy  is  now  very  simply 


P  >  P  +  P 
total  -  refl  rad 


(24) 


with  the  eqmlity  sign  holding  if  and  only  if  we  are  discussing  a  lossless  structure. 
For  at  least  part  of  otir  discussion,  we  shall  not  make  any  such  restriction;  our 
equations  above  are  perfectly  valid  even  if  absorption  is  taking  place  within  the 
structiire,  so  long  as  the  system  remains  linear  in  the  general  sense  (in  the  sense 
of  superposition  holding  for  all  voltages  and  currents).  We  may  also  note  that 
althoiigh  for  any  linear  passive  network,  the  S  matrix  will  be  symmetric  (S.  .  ■  S  .  ) 
by  the  le.w  of  reciprocity,  our  forms  above  have  not  Included  this  requirement,  and 
our  equations  include  the  possibility  of  non-reciprocal  systems  (non-syametric  S), 
provided,  again,  they  are  still  linear  in  the  sense  of  superposition  remaining  valid. 
With  these  comments  in  mind,  we  finally  rewrite  Eq.  (24)  In  terms  of  Eqs.  (l8), 

(19),  and  (23),  as* 

x'^x  >  x^S^Sx  +  x"^  p  X  (25) 

and  this  equation  must  hold  for  amy  arbitrary  choice  of  the  excitation  voltages,  x. 


* 

I  am  indebted  to  my  colleague.  Dr.  J.  E.  Storer,  for  pointing  out  the  essential 
matrix  nature  of  the  fundamental  eqxiations,  following  upon  my  original  derivations. 
Many  of  the  results  below  would  not  have  been  obtained  without  this  recognition, 
and  the  resulting  compactness  of  the  mathematics.  The  general  fundamental  result 
in  Eq.  (39)  was  pointed  out  by  Dr.  Storer,  immediately  upon  this  recognition. 

He  also  pointed  out  the  possibility  for  extending  the  studies  to  discussions  of 
the  scattering  matrix,  thus  stimulating  and  contributing  to  the  further  investiga¬ 
tion  reported  in  the  canpanion  paper,  Psurt  II.'  ' 


-9- 


ARM-237 


We  may  note  that,  if  we  introduce  the  identity  matrix,  I, 

we  can  rewrite  Eq.  (25)  as 

x"'(l  -  S'^S  -  r)  x>  0  (27) 

end,  since  this  holds  for  arbitrary  x,  it  states  that  (l-S^S-  P)  must  be  a  posi¬ 
tive  semi -definite  matrix  for  a  lossy  network,  and  must  be  the  zero  matrix  if  and 
only  if  a  lossless  struct\ire  is  involved.  As  will  be  shown  Immediately,  it  is  also 
Hermitian,  and  this  is  sufficient  to  provide  much  valmble  information.  However, 
we  will  ultimately  wish  to  return  to  Eq.  (25),  and  hence  will  continue  to  use  it 
directly  throughout  our  derivations. 


5.  Limitation  on  Radiation  Efficiency  Implied  by  the  Conservation  of 
Energy  Equations; 

Let  us  first  note  that  both  P  and  the  matrix 
W  =  S'^S 


are  Hermitian,  which  is  to  say  that 


r+  =  r 


^kj  “  ^  kj 


\t  =\l 


The  second  of  these  follows  directly  from  Eq.  (28), 


-  (S'^S)'^  =  (S'^){S'^)'^  =  S'^S  =  W 


(28) 


(29) 
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euid  the  first  from  the  definition  of  T  in  Eq.  (22),  and  the  symmetry  prqperty  of  P 
In  Eq.  (9), 


/  +\  * 

^  ^kj  “  ^Jk 


^Jk  ^k)* 


^Jk  \ 


\  "kj  ‘Ij 


kJ 


Secondly,  both  f  and  W  are  positive  semi-definite,  i.e.,  for  any  vector  z, 
z"*"  r  z  >  0 

.  (30) 

z  Wz  >  0 

The  second  of  these  follows  almost  Immediately  again  from  Eq.  (28) 
z^Wz  =  z^S^Sz  =  (Sz)^(Sz)  >  0 


+  M 

fi.e.,  for  any  colimin  vector  t,  t  t  =  Y  t  >  0,  with  equality  only  if  all  t.  =  O). 

k=l  ' 

The  first  follows  likewise  by  noting  from  Eqs.  (22)  euid  (8), 

z  rz-  I  . 

k,J=l 


M 


I,  \  \(»> 


dO  >  0 


since  the  product  in  the  Integral  is  now  the  square  of  the  absolute  value  of  the 
magnitude  of  a  con5)lex  vector.  Now,  with  both  W  and  P  known  to  be  positive  seml- 
definlte  and  Hermltian,  we  know  from  matrix  theory  the  inqportant  property  that 
each  has  exactly  M  positive  real  or  zero  eigenvalues,  and  M  corresponding  mutually 
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It  is  in^jortant  to  note  that  each  eigenvector  of  a  matrix  is  determined  only  to 
within  a  scale  factor.  The  normalization  to  unit  length,  i.e.,  the  requirement 

z'*^z  «  1,  is  necessary  if  the  eigenvectors  are  to  he  used  to  construct  a  unitary 
matrix  with  the  properties  given  in  Eqs.  (31)  and  (32).  However,  the  coa^ionents 
of  z  are  still  conqjletely  Indeterminate  to  within  a  common  phase  factor  of  the 

itfj 

form  e  .  Stated  another  way,  suppose  we  consider  any  matrix  D  which  is  diagonal, 
with  elements. 


k  A  J 


i.e.,  D  of  the  form 


where  {6^,  are  arbitrary  real  numbers.  Then  it  is  reeidlly  verified  that 

any  such  D  itself  is  \mitary  and  that  any  matrix  V  of  the  form 
V  =  UD 


is  also  unitary,  and  also  diagonalizes  F  into  7.  Forms  like  D  are  encountered  in 
Part  II  in  attempting  to  characterize  the  form  of  the  scattering  matrix. 
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If  we  return  to  Eq.  (25),  and  introduce  U  as  Just  defined  by  transforming  from 
X  to  a  new  (equally  arbitrary)  vector  z, 

X  =  Uz 


+ 

X 


+Tt+ 

z  U 


(34) 


we  ceua  rewrite  Eq.  (25)  as 

z'^U'^Uz  >  z'^U'^S'^SUz  +  z'^U'*’  r  Uz 
Or,  using  Eqs.  (31)  and  (32) 

z^z  >  z^U^ S^Uz  +  z^ T  z  (35) 

We  may  rewrite  this  inequality  as 

z^(l-‘?')z  >  z'^Tz  (36) 

where 

T  -  uVsU  (37) 

is  readily  observed  to  also  be  Hermitian,  and  positive  semi-definite.  Now  tbe 
right-hand  side  of  Eq.  (36)  is  positive  or  zero,  for  arbitrary  z,  while  on  the 
left-hand  side,  the  matrix  (l-  bas  non- zero  elements  only  on  the  diagonal.  It 
follows  frcsn  the  arbitrariness  of  z  under  which  Eq.  (36)  holds,  that  these  diagonal 
elements  of  (l-  ^)  must  hence  individually  be  positive  or  zero,  i.e.,  we  must  be 
able  to  write 


(38a) 

4  -  ^  -  ’'k  i  ° 

(38b) 
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This  is  an  important  resiilt.  In  particular,  Eq.  (38b)  implies  that  the  largest 


of  the  eigenvalues  of  r  cannot  exceed  unity, 

SI  (35) 

If  we  recall  the  form  of  i.e., 

^  kj  “  \  “Ij 

where  the  P.  .  are  fixed  by  the  design  of  the  beams,  we  see  that  Eq.  (39)  places 

direct  1 imitatlnnH  on  the  values  of  the  beam  radiation  efficiencies,  Iq^l^.  In  a 

(1) 

con5)anion  paper,  we  derive  further  from  the  conservation  of  energy  equations, 
some  of  the  properties  of  the  scattering  matrix,  S.  In  the  remainder  of  this  paper 
we  will  discuss,  via  exan5)les,  the  physical  meemingfulness  of  the  result  in  Eq.  (39)- 

6,  Example;  Equal  Beam  Efficiencies  (Angle  Diversity  Receiving  Arrays) 

An  important  and  sinqile  result  is  obtained  from  Eq.  (39)  for  the  case  when  all 
the  beams  are  designed  to  be  a  priori  equal,  i.e. 

qj^  •=  q  for  all  k  (^) 

This  has  been  the  usual  design  case,  for  example,  in  angle  diversity  receiving  arrays 

of  the  kind  under  recent  study  in  connection  with  tropospheric  scatter  reception.; 

'  ’  A.B.  Crawford,  D.C.  Hogg,  and  W.H.  Kummer,  "Studies  in  Tropospheric  Propagation 
Beyond  the  Horizon,"  BSTJ,  IO67-IIT8  (1959). 

(o') 

J.H.  Vogelman,  J.L.  Ryerson,  and  M.H.  Bickelhaupt,  "Tropospheric  Scatter  System 
Using  Angle  Diversity,"  Proc.  IRE,  pp.  688-696  (1959)* 

R.  Bolgiano,  Jr.,  N.H.  Bryant,  and  W.E.  Gordon,  "Diversity  Reception  in  Scatter 
Communications  with  Emphasis  on  Angle  Diversity,"  Cornell  Univ.,  Final  Report, 
Part  I,  on  AF30(602)-1717,  Jan.,  I958. 

C.  Chu,  A.H.  Wren,  and  J.  LaRue,  "Evaluation  of  the  Pincushion  System,"  Univ. 
Mich.  Rad.  Lab.  Kept.  2872-l-T,  RADC-TN-6O-5O,  Feb.  I96O. 

W.R.  Richard  and  M.H.  Bickelhaupt,  "Multiple  Angle  Diversity  Design  Considera¬ 
tions,"  Rome  Air  Development  Center,  RADC-TN-60-22,  April,  196O. 
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and  it  was,  in  fact,  difficulties  which  were  first  noticed  d\aring  a  related  study 

which  originally  motivated  the  presait  investigation.  In  disucssing  a  receiving 

eirray,  rather  than  a  transmitting  array,  we  are  invoking  the  well-known  reciprocity 

relationship  between  the  receiving  and  transmitting  antennas,  which  can  be  readily 

shown  to  apply  to  the  radiation  efficiency  defined  by  above.  I.e.,  if  all  M 

beams  are  allowed  to  receive  simultaneously  with  matched  receivers,  one  may  calculate 

the  actual  power  received  by  each  receiver  by  first  calculating  the  power  which  would 

be  received  according  to  its  corresponding  beam  pattern,  as  if  no  other  beams  existed, 

I  1 2  'll' 

and  then  multiplying  by  the  radiation  efficiency  for  that  beam. 

Returning  to  Eq.  (Uo),  we  see  that  it  implies 


(41) 


It  theii  reeuiily  follows  that  the  ,  the  eigenvalues  of  P ,  are  directly  related 

to  ,  the  set  of  eigenvalues  of  P,  by 

5,,  (42) 


(Y) 

'  'D.E.  Johansen  emd  S.  Stein,  "Theory  of  Antenna  Performance  in  Scat+er-Type 
Reception."  to  be  published  in  IRE  Trans.  PGAPj  also  see  report-  under  same 
title  by  Stein,  Johansen,  and  Starr,  Hermes  Electronics  Co.,  AP.CRC-TR-59-191^ 

30  September  1959* 

A  note  of  warning:  In  many  of  our  results,  there  is  no  limitation  on 
whether  or  not  the  eintenna  system  is  reciprocal.  However,  whenever  we  discuss 
a  receiving  array  in  terms  of  transmitting  properties,  the  entire  antenna  system 
is  being  assumed  to  be  reciprocal.  In  an  actual  case  of  a  non-reciprocal 
antenna  system,  one  will  generally  want  to  consider  the  receiving  properties  in 
terms  of  an  equivalent  reciprocal  transmitting  antenna}  where,  however,  it  may  then 
no  longer  be  possible  to  relate  receiving  properties  to  measurements  made  with 
the  real  antenna  system  transmitting. 
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and  Eq.  (39)  becomes 


(^3) 


k^majc 

IbuB,  in  this  case,  the  inherent  limitations  due  to  conservation  of  energy  are 
conpletely  related  to  the  array  of  beam  cross-coupling  factors.  Incidentally,  it  is 
well  known  that  for  any  Hermitian  matrix,  the  sum  of  the  eigenvalues  equals  the  sum 
of  the  diagonal  elements  (the  trace  of  the  matrix).  Th\xs  the  largest  eigenvalue 
exceeds  the  average  of  the  diagonal  elements,  so  that  since  all  diagonal  elements  of 


P  have  the  same  value  *  1,  |q|  <  1  is  iit5)lied;  and  |q|  =  1  is  possible  only 

if  adl  eigenvalues  are  equal,  which  in  t\irn  is  readily  observed  to  require  that  all 
non-diagonal  elements  of  p  vanish,  i.e.,  no  beams  overlap  (e.g.,  see  the  next  sec¬ 
tion  ) . 


In  Appendix  1,  we  discuss  some  specific  numerical  examples  which  indicate  the 
order  of  magnitude  of  the  Iosl  in  radiation  efficiency,  for  what  we  regard  as  some 
very  realistic  antenna  beam  patterns  for  angle  diversity.  In  particular,  we  note 
that  although  previous  angle  diversity  analyses  have  considered  the  effects  of  beam 
overlap  with  respect  to  correlation  in  observed  propagation -induced  fluctuations  at 
the  beam  output  receivers  (see  also  Section  8),  they  have  generally  ignored  the 
(basic8j.ly  independent)  question  of  the  \mavoidable  "insertion  losses"  in  the  antenna 
Junction  whose  existence  has  been  shown  above.  A  general  statement  Justifying  such 
an  approach  seems  to  be  the  assertion  that  if  beams  are  squinted  so  as  to  cross  over 
at  the  3  db  points  or  lower,  then  in  effect  they  represent  essentially  independent 
receptors  of  power.  The  inconsistency  of  this  assiimption  is,  we  feel,  evident  frcsn 
Appendix  1;  the  numerical  results,  nevertheless,  we  found  to  be  somewhat  startling. 

In  addition,  we  refer  the  reader  interested  in  angle  diversity  to  the  further 
discussion  in  Section  8,  cmd  the  interesting  numerical  exoo^le  and  discussion  in 
Appendix  3- 
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7.  Example;  Two  Limiting  Cases 

As  further  relevant  exanrples,  we  consider  two  special  cases: 
(a)  No  beam  overlap. 

In  this  case 


J  It 


Thus,  for  this  case 


k  -  J 


0  k  ^  J 


(44b) 


and,  by  inspection,  the  eigenvalues  of  P  are  exactly 


(44c) 


Th\is  if  all  q 


( |q 


are  equal,  all  beams  can  be  made  to  have  100%  efficiency 
i),  as  one  might  expect  physically. 


A  further  interesting  point  now  arises  from  this  example.  Suppose  even  though 
the  beams  are  unco^qjled,  the  design  aims  at  different  ^or  each  beam.  Then  one 

or  more  beams  corresponding  to  the  largest  value  will  have  100^  efficiency,  but  the 
others  will  not.  That  is,  although  potentially  all  beams  could  be  made  100^  efficient 
in  this  non-overlap  situation,  a  poor  a  priori  choice  of  the  relative  qj^  would  iii5>ly 
building  into  the  system  cross-couplings  which  are  not  "necessary".  This  may  appear 
to  be  a  very  peculiar  kind  of  statement;  nevertheless,  it  appears  that  there  are 
many  more  general  antenna  system  engineering  problems  where  such  a  priori  choices 
must  be  made  (usually  the  choice  being  that  all  qj^  sure  eqvial,  perhaps  without 
realization  that  alternatives  exist).  A  further  discussion  of  this  issue  is  given 


in  the  next  section. 
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In  addition,  in  Appendix  2,  we  indicate  how  one  might  extend  these  non-overlap 
results,  via  perturbation  theory,  to  a  case  of  much  practical  design  interest,  namely 
that  in  which  the  amount  of  beam  coupling  between  all  beams  is  small. 

(b)  Coitplete  Beam  Overlap 

From  a  mathematical  point  of  view,  it  is  illviminating  to  consider  the 
other  extreme  limiting  case  where  all  beams  overlap  completely.  Physically,  this 
corresponds  to  having  several  "taps"  at  some  point  down  the  line  from  a  single  feed. 
That  is,  we  really  do  not  have  separate  beeuns,  and  certainly  cannot  extract  at  each 
tap  an  amount  of  power  equal  to  the  total  power  received  by  the  beam;  yet  we  insist 
that  a  correct  mathematical  formulation,  such  as  we  claim  to  now  have,  must  be 
capable  of  fully  and  correctly  accounting  for  this  situation  (this  sinple  limiting 
exan^jle,  again,  being  one  of  the  sotirce  points  for  the  entire  investigation). 

For  complete  overlap,  we  have 

for  all  j,  k  (14.5) 

It  is  readily  determined  that  the  eigenvalues  of  P  conqjrise  an  (M-l)-fold  degenerate 
zero  eigenvalue,  and  one  with  the  value 

r  .  X 

J 

Th\is,  without  discussing  the  individual  beam  radiation  efficiencies,  we  see  immedi¬ 
ately  from  Eq.  (39)  that  we  have 
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which  is  just  what  o\ir  physical  intuition  requires.  While  there  is  nothing  really 
startling  in  our  achieving  this  result,  since  we  ass\ime  o\ir  mathematics  have  been 
correct,  it  is  quite  satisfying  to  know  that  the  formulation  routinely  describes 
this  limiting  case,  as  well  as  other  cases  not  quite  so  simple.  We  are  not  aware 
of  any  previous  mathematical  description  for  this  general  problem  which  has  this 
satisfying  property  between  all  extremes. 


8.  Exampl e i  Specification  of  Relative  Radiation  Efficiencies  (Application  to 
Optlmm  Angle  Diversity) 


Let  us  return  to  the  general  result  of  Eq.  (39)>  and  consider,  for  a  general 
form  of  the  matrix,  P,  the  question  of  specifying  relative  radiation  efficiencies. 

To  emphasize  what  is  being  specified,  it  is  illuminating  to  introduce  these  relative 
values  as  dimensionless  ratios,  K^,  and  to  assume  that  all  the  q^  are  now  in  the 
form 


=  \<i 

The  elements  of  r  are  thus  in  the  form 


(45a) 


'•ij  =  bf  Vi/j 

-  bf  Hu  <‘‘5'=) 

Thus,  if  we  can  find  the  eigenvalues,  oi  the  matrix  ^ ,  we  will  have 

Y,  =  |ql=  Cea) 

and  from  Eq.  (39) 


k  max 


(46b) 
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At  this  time,  this  much  is  straightforward;  however,  we  regard  the  availability 
of  an  arbitrary  pre-specification  of  the  as  constituting  an  extremely  interest¬ 
ing  problem  area.  It  appears,  moreover,  that  it  is  precisely  with  this  problem 
that  one  necessarily  makes  contact  with  the  question  of  the  use  of  the  multiple- 
beam  antenna  system.  In  many  ajqilications ,  such  as  angle  diversity  referred  to 
above,  or  the  case  of  a  phased  array  in  which  only  one  beam  is  active  at  a  time, 
the  question  of  the  radiation  efficiency  of  a  single  beam  and  of  Individual  line 
cross-couplings  is  i)aramo'unt.  In  other  applications,  more  than  one  transmitting  beam 
may  be  excited  at  a  time,  and  it  is  then  the  performance  of  the  superposition  of  coupled 
currents  which  is  relevant.  Further,  for  example,  in  examining  angle  diversity  for 
tropospheric  scatter  systems,  one  would  want  to  choose  the  to  maximize  the  diver¬ 
sity  performance;  in  other  systems,  one  might  want  to  impose  a  constraint  of  siaxi- 
mum  (or  minimum)  antenna  system  directionality  in  certain  directions;  in  others, 
an  approximation  within  specified  tolerances  to  some  conposite  beam  i)attem,  with 
maximum  power  transmitted  or  received,  ml^t  be  the  goal. 

That  this  problem  of  the  optimum  choice  of  the  is  Important  we  have  already 
seen  from  our  discussion  of  the  zero  coupling  case  in  the  previous  section.  It 
was  pointed  out  that,  if  all  beams  are  decoupled,  smy  choice  of  unequal  would 
lead  to  less  than  100^  radiation  efficiency  for  some  of  the  beams,  but  that  in 
this  fully  decoupled  case  there  was  no  inherent  reason  for  making  such  a  poor 
choice.  While  for  this  particular  problem  the  correct  choice  was  obvious  on  phy¬ 
sical  grovinds  (see  also  the  discussion  of  the  weak  coupling  case  in  Appendix  2), 
it  is  not  at  all  apparent  as  to  what  might  be  lost  in  a  more  general  situation 
where  an  arbitrary  choice  (say  all  equal)  is  made  for  the  K^.  Thus  the  general 
unsolved  problem  is  to  learn  more  about  the  implications  of  the  relative  values  of 
the  q^,  as  they  affect  the  Y^,  and  in  turn  through  the  latter  (the  conservation 
of  energy  constraints),  the  absolute  scale  on  all  the  q^. 
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As  an  interesting  exan^jle  of  vhat  such  a  study  might  lead  to  in  a  peu:ticular 
configuration,  we  can  formulate  the  problem  in  the  angle  diversity  case.  (Unfor- 
txuiately,  we  cannot  solve  the  problem  in  general.  Nevertheless,  as  we  indicate 
later  by  an  example  in  Appendix  3>  in  many  particular  ernes,  solution  may  be  pos¬ 
sible  without  exorbitant  effort.)  The  quantity  of  direct  Interest  in  these  pro¬ 
blems  is  either  the  probability  distribution  of  the  fluctuations  in  received  power 
or  some  closely  related  system  performance  statistic.  The  basis  of  angle  diver¬ 
sity  as  applied  to  the  fading  signals  of  tropospheric  scatter  is  the  desire  to  talce 
advantage  of  experimental  observations  that  the  instantaneous  angle  of  arrival  of 
the  waves  appears  to  wander.  A  number  of  receiving  beams  are  thus  used,  each  point¬ 
ing  in  a  different  direction,  emd  their  outputs  combined  via  an  appropriate  network. 
The  best  combining  rule,  if  thermal  noise  fluctuations  are  independent  from  receiver 
to  receiver,  is  so-called  maximal-ratio  combining,  the  effect  of  which  is,  assuming 
equal  noise  in  all  receivers,  to  edd  the  effective  powers  of  all  the  receivers. 

The  output  noise  level,  on  the  a88uii5)tlon  made,  is  the  same  as  that  in  euay  single 
receiving  channel.  Further,  assuming  that  the  fluctuations  in  signal  are  described 
by  so-called  Rayleigh  envelope  fading,  possibly  correlated  between  receivers  due  to 
beam  overlap,  the  characteristic  function  (Laplace  transform  of  the  probability 

/Q) 

density)  has  the  form'  ' 

Here  z  is  the  transform  variable,  while  L  is  the  M  x  M  conplex  covariance  matrix 
among  the  cemplex  characterization  (a  band-pass  system  is  assumed)  of  the  fluctiiation 

^®^J.N.  Pierce  and  S.  Stein,  "Multiple  Diversity  with  Nonindependent  EBidlng," 

Proc.  IRE,  48,  pp.  89-104  (i960). 
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process  at  each  receiver.  In  the  paurlaace  of  our  present  investigation. 


where  the  bar  indicates  a  statistical  average  over  the  fading  fluctuations,  and 

each  Bj  is  the  (ccai^jlex)  voltage  (effective  rms  into  a  1  ohm  load)  which  would  be 

til 

obtained  at  the  receiver  output  if  the  j  beam  were  received  with  its  specified 
pattern,  with  no  other  beams  present  in  the  configuration  (no  cross -coupling  to 
other  beams ) .  It  has  further  been  pointed  out'  '  that  for  non-coherent  FSK  digital 
data  communications,  ^  F(z)  in  Eq.  kj  gives  the  average  bit  error  rate,  if  z  is 
taken  to  be  z  =  ,  where  is  the  noise  level.  Better  yet,  if  all  signal 

levels  are  normalized  to,  say,  a  median  level,  z  can  be  taken  as  proportional  to 
the  median- signal,  noise  ratio  (see  Appendix  3)*  Hence  F(z)  in  Bq.  (47)  may  in 
itself,  for  a  specific  z,  be  considered  a  useful  performance  measure  for  the  sys¬ 
tem.  In  writing  Eqs.  (47)  and  (48),  there  is  no  presunqption  that  the  average  level 
of  the  fading  signals  in  the  respective  receivers  (the  e^)  are  the  same;  in  general, 
depending  upon  the  nature  of  the  transmitting  antenna  pattern  and  the  actual  receiv¬ 
ing  array  configuration,  these  average  levels  will  differ.  Finally,  under  very 

(7) 

reasonable  assim^jtions' '  ,  it  is  permissible  to  assume  that  one  can  define  an 
apparent  distribution  of  fluctuating  sources  "in  the  sky"  equivalent,  in  reception, 
to  the  actual  fields  set  up  by  the  transmitting  antenna.  In  terms  of  the  intensity 
distribution  (mean  square  value  at  each  direction)  of  this  equivalent  source  dis¬ 
tribution,  one  can  write 


^  j  R^(0)  •  Rj(0)  g(0)  dO 


where  g(o)  is  the  aforementioned  intensity  distribution  and  the  are  the 

normalized  beam  patterns  defined  earlier  by  Eqs.  (4)  and  (5).  If  we  thus  define 
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an  array  of  "beam  correlation  factors"  by  the  equation 

j  *^®  »  ^50) 

we  can  write  the  matrix  L  in  the  form 


~  ®kj 


(51) 


(quite  analogous  to  the  form  of  the  r  -matrix), 
dimensionless  (Eq.  45a),  and  write 

'  hr  <  %  ■'j 


We  can  further  introduce  again  the 


(52) 


Thus  the  over-all  optimization  of  the  angle-diversity  system,  assiming  the  treuis- 
mitted  field  and  the  receiving  beams  to  be  alreeuiy  specified  and  fixed,  consists 
in  finding  that  set  *  K  \  which,  along  with  the  associated  best  possible  value  for 
jqj  (Eqs.  45-46),  minimizes  F(z)  in  Eq.  (47)  for  some  specified  z.  The  implica¬ 
tion  of  this  statsnent  of  the  problem  is  that  with  a  specified  transmitted  field 
and  array  of  receiving  beams,  it  might  be  optimum  to  enhance  the  efficiencies  of 
some  of  the  beams  at  the  cost  of  degrading  others,  to  achieve  an  over-all  optimiun 
performance.  The  formulation  above  seems  to  provide  a  precise  mathematical  state¬ 
ment  of  a  problem  which  has  been  remarked  upon  only  superficially,  at  best,  in 
previous  studies.  The  solution  to  this  problem,  with  further  optimization  of  the 
receiving  beam  array  configuration  and  nature  of  the  transmitted  field,  seems 
likely  to  provide  useful  information  as  to  the  true  effectiveness  of  angle  diver¬ 
sity  systems. 

In  Appendix  3>  ve  consider  a  typical  angle  diversity  receiving  array  by  way 
of  illustration,  with  some  Interesting  conclusions  (at  least  for  the  numbers  chosen). 
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9.  Conclusions  and  Unsolved  Problems 

A  general  formulation  has  been  given  vhich  points  up  the  in5>ortance  suid  some 
of  the  iii5)lications  of  conservation- of- energy  relationships  in  mviltiple-beam  antenna 
systems.  These  inpli cations  appear  as  limitations  on  radiation  efficiencies  for 
individual  beams.  In  Part  II, further  in^ilicatlons  are  given  in  terms  of  char¬ 
acterizing  the  form  of  the  scattering  matrix,  i.e.,  the  eurray  of  feed-line  self- 
and  cross -coupling  factors.  These  results,  together,  appean  quite  importsmt  for 
elucidating  inherent  antenna  system  design  limitations  following  upon  particular 
choices  of  beam  patterns  in  such  multiple-beam  systems. 

It  appears  that  further  useful  results  remain  to  be  obtained.  These  unsolved 
problems  appear  to  lie  in  better  understanding  of  certain  matrix  problems,  and  In 
applying  maximum  ingenuity  towards  analyzing  particular  classes  of  matrices.  With 
respect  to  the  material  discussed  here  in  Part  I,  the  major  aspect  which  we  regard 
as  unsolved  is  the  significance  of  the  a  priori  choices  of  the  relative  beam  radia¬ 
tion  efficiencies  (e.g.,  the  defined  in  Section  8).  In  particular,  the  open 
question  is  whether  one  can  determine  optimum  choices,  using  relevant  criteria  as 
to  the  desired  system  performance.  Even  in  the  angle  diversity  case  cited,  where 
the  problem  has  been  reduced  to  a  specific  mathematical  statement,  one  does  not 
have  a  general  solution.  For  other  systems,  the  deduction  of  the  mathematical 
statement  from  vhich  c^tlmizatlon  can  be  determined  may  represent,  in  each  case, 
a  fomldable  problem  in  itself. 
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GAUSSIAN  BEAM  PATTERNS 
TAPERED  lUUMINATION  p:| 
TAPERED  ILLUMINATION  p:2 
UNIFORMLY  ILLUMINATED  APERTURE 
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9;  HALF-SQUINT  ANGLE  BETWEEN  BEAMS 


Figure  3.  Beam-Coupling  Factors 
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Figure  4.  Radiatr 
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Figure  7.  Error  Rate  Versus  |  K 
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APPENDIX  1;  A  NUMERICAL  EXAMPLE  FOR  RADIATION  EFFICIENCY* 

In  Figure  2,  we  show  three  beam  patterns  of  interest  (field  patterns)  for  a 
circular  aperture  with  linear  polarization.  Two  of  these  are  idealized  patterns: 
These  are  for  the  "\aniformly  Illuminated"  aperture  of  radius  a,  and  the  gaussian- 
illuminated  aperture  (with  the  aperture  currents  extending  to  infinity).  The 

third  is  a  tapered-aperture  distribution  of  the  form  (1 - ^  J  ,  with  p  =  1.  For 

all  three,  the  parameters  of  the  antenna  were  chosen  so  that  the  main  lobes  are 
essentially  of  identical  width  at  the  3  db  points.  The  tapered  distribution  r^re- 
sents  a  good  approximation  to  actual  designs  which  aim  at  lower  side-lobe  levels 
than  associated  with  the  imiformly-illumlnated  aperture,  while  at  the  same  time 
avoiding  the  losses  due  to  spillover  (which  do  not  appear  in  these  theoretical 
calcvilations  ) . 

We  assume  next  that  for  each  type  of  beam,  a  multiple -beam  system  is  con¬ 
structed  with  a  common  aperture,  with  the  various  "squint"  angles  (directions  of 
pointing)  obtained  by  appropriate  uniform  phasing  of  the  currents  across  the  aper¬ 
tures.  With  this  assunption,  we  can  readily  compute  the  beam-coupling  factors 
between  any  pair  of  beams,  (incidentally,  the  calculations  are  most  readily  accom¬ 
plished  by  recalling  the  Fourier  Transform  relationship  between  beam  patterns  and 
corresponding  aperture  distributions  j  and  then  transforming  Eq.  (8)  as  a  special 
case  of  a  convolution  integral  to  give  an  equivalent  integral  over  the  aperture 
plane.)  In  Figure  3,  we  have  plotted  as  a  function  of  the  squint  single  between  suiy 

pair  of  beams  of  like  type,  the  associated  beam-coupling  factor,  P.  An  additional 
_ 

I  am  indebted  to  my  colleague,  Mr.  W.  J.  Cowan,  for  providing  me  with  the 
material  of  this  section,  especially  Figures  2,  3,  and  4;  and  as  well,  the 
numerical  data  used  in  the  example  of  Appendix  3* 
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curve  has  been  added  here  for  the  p  «=  2  taper,  whose  beam  pattern  would  be  even 
closer  to  the  gauss Ian  In  Figure  2  than  the  p  =  1  pattern.  The  abscissae  in 
Figure  3  are  half  the  total  squint  angle  between  beams,  with  the  same  angular  scale 
as  in  Figure  2.  Thus,  for  beams  crossing  over  at  the  3  db  points,  we  see  by  compar¬ 
ing  Figures  2  and  3  that  the  overlap  factor  varies  from  0.14-5  for  gausslan  shaped 
beams,  to  about  0.35  Tor  the  tapered  illuminations,  to  about  0.1  for  the  uniformly 
illuminated  case. 

For  two  beams,  with  eqvial  a  priori  radiation  efficiencies,  the  eigenvalues  of 
interest  are  those  of  the  P-matrix,  obtained  by  solving 


Or 


1  -  X  P 

P  1  -  X 


=  1  -  P 


Xg  =  1  +  P 


(A1.-1) 


From  Eq.  (43),  the  associated  upper-bound  radiation  efficiency  is  thus  given  as 


max 


1  +  P 


(AI-2) 


For  P  =  0.49, 
P  =  0.3  , 
P  -  0.14, 


q  =  0.67 

^  max  ' 


max 

12 

max 


0.76 

0.88 


(Al-3) 


These  are  shown  as  the  case  M  =  2  in  Figure  4.  It  is  interesting  to  note  that 
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with  only  the  two  beams,  t}ie  efficiencies  are  already  reduced  by  25^,  or  about 
1.4  dbl  In  Figure  4,  we  have  also  Indicated  additional  radiation  efficiencies 
for  particular  higher-order  configurations  for  which  calculations  have  been  made 
for  other  purposes.  In  all  cases,  equal  a  priori  radiation  efficiencies  are  assumed 
(also  see,  however.  Appendix  3)*  The  circles  in  the  sketches  of  the  configurations 
show  how  the  -3  db  contours  intersect  for  the  assumed  configurations.  It  is  inter¬ 
esting  to  note  that  the  data  do  not  seem  to  lie  on  smooth  curves.  This  is  parti¬ 
cularly  true  for  the  M  ~  5  case.  We  attribute  this,  comparing  the  sketches  of  the 
configurations,  to  the  fact  that  it  is  the  actual  couplings  of  each  beam  to  adjacent 
beams  which  is  Inportant,  and  these  are  especially  different  for  the  M  =  5  configu¬ 
ration  than  for  either  the  M  =  7  or  M  =  l8  cases. 

Incidentally,  the  resiilts  in  Figure  3  indicate  quite  different  beam  overlap 
factors,  and  hence  quite  different  cross-coiplings,  for  beams  which  differ  only  in 
their  side  lobe  structure  but  have  almost  identical  main  lobes.  In  one  sense,  this 
points  out  the  importance  of  the  side-lobe  structure  in  multiple -beam  systems.  In 
another  sense,  it  verifies  a  common  engineering  fact,  that  one  may  make  minor  changes 
in  a  feed  structure  which  essentially  leave  the  main  beam  xmchanged  (presumably, 
however,  radically  altering  the  side-lobes),  while  radically  reducing  cross -coupling 
effects.  It  is  also  interesting  to  note  from  Figure  3  that  those  beam  patterns 
(tapered  illuminations)  which  more  closely  approach  the  ideal  of  a  monotonic  decrease 
of  field  off  the  direction  of  pointing,  are  also  those  which  are  worst  in  terms  of 
producing  beam  overlap  and  hence  cross -coupling  effects.  In  any  case,  beam  overlap 
effects  appear  to  be  far  from  insignificant. 
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APPENDIX  2:  WEAK  BEAM  OVERLAP  -  A  PERTURBATION  CALCULATION 


In  this  case,  ve  assume  for  the  overlap  factors. 


kj 


«  1  ,  J  k 


(A2-1) 


An  array  of  overlap  factors  of  this  type  is  extremely  in^jortant  from  an  engineering 
viewpoint,  since  it  is  likely  to  represent  the  situation  where  some  overlap  is 
occurring,  but  every  atten^jt  has  been  made  in  design  planning  and  because  of  opera¬ 
tional  requirements  to  keep  it  small.  For  particular  cases,  where  conditions  of 
symmetry  and  the  like  sufficiently  simplify  the  problem,  an  exact  solution  may  be 
possible.  This  is  shown,  for  exainple,  by  the  illustrative  problem  treated  in 
Appendix  3.  In  this  section,  however,  we  wish  to  indicate  a  more  generally  appli¬ 
cable  method  for  obtaining  an  approximate  solution  in  a  general  case,  by  applying 
standard  perturbation  techniques,^^^^^^^^  with  the  zero-coupling  case  in  Eq.  (44) 
regarded  as  the  unperturbed  case.  I.e.,  we  write  the  P  -matrix  for  the  weak- 
coupling  situation  in  the  form 


r  =  +  €  r’  (A2-2) 

where  p  is  the  matrix  for  the  zero-coupling  case,  p'  is  a  matrix  whose  ele¬ 
ments  are  of  the  same  order  as  those  of  ajid  e  is  a  perturbation  parameter  for 

which  it  is  assumed  that 

€  «  1  (A2-3) 

L.I.  Schlff,  Quantum  Mechanics,  McGi^w-Hill,  New  York,  1949,  p.  149-156. 

Bellman,  Introduction  to  Matrix  Analysis,  McGraw-Hill,  New  York,  I96O,  p.  61-63. 
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Thus,  referring  to  Eq.  (^4-4),  we  have 


i(0)  ^ 


=  6 


Jk  I'^k 


(A2-4) 


Further,  for  we  can  formalize  the  results  in  Eq.  (44),  by  describing  its 


th 


eigenvector,  and  corresponding  m  eigenvalue,  7^,  as  satisfying  the  matrix 
equation 


,(o) 


u  =  7  u 
m  mm 


(A2-5a) 


or,  in  ccagjonent  form, 

^  =  7„u^  (A2-51.) 

where  u^^  is  the  k  conqjonent  of  the  m  "  eigenvector.  Referring  to  Eq.  (44),  we 
readily  observe  that 


7^  =  |q^|^  (A2.6a) 

and  that  the  detailed  con5>onents  of  the  corresponding  normalized  m  ”  eigenvector  have 
the  very  simple  form 


u 


mk 


0  ,  k  / 


m 


u  =1 
mm 


(A2-6b) 


In  addition  to  all  eigenvectors  being  normalized  to  length  unity,  we  may  recall  that 
for  a  Hermitian  matrix,  they  also  form  a  mutually  orthogonal  set.  For  the  case  of  a 
degenerate  eigenvalue,  there  is  some  arbitrariness  in  the  construction  of  this  mutually 
orthogonal  setj  we  will  discuss  this  point  further  below  in  connection  with  the  per¬ 
turbation  calculation  (we  have  seen  already  that  we  may  indeed  have  much  physical 
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interest  In  a  case  vhere  Boine  or  aU  the  are  equal  In  magnitude,  in  vhich  event 
the  M  X  M  matrix  has  vrp  to  an  M-fold  degeneracy). 

The  perturbation  calculation,  in  effect,  seeks  to  observe  the  (presumed)  gradual 
ajid  continuous  changes  in  the  eigenvalues  and  eigenvectors  of  T  as  the  perturbation 
parameter,  €,  in  Eq.  (A2-2)  increases  slowly  from  zero.  Since  the  diagonal  elements 
of  any  P  -matrix  are  always  Just 


p(0) 

’  kk 


(A2-7) 


we  see  that  any  perturbing  matrix,  p’,  in  Eq.  (A2-2)  must  have  identically  zero 
elements  along  the  main  diagonal. 


Pj^  =  0  for  all  k  (A2-8a) 


Furthermore,  although  €  is  to  be  a  parameter  in  the  perturbation  calculation,  we 

will  always  wish  to  choose  as  oxir  point  of  calcixlation  the  combination  of  e  and  p’ 

(0) 

such  that  for  the  off-diagonal  elements  (recalling  that  Pj^^  =  0,  J  /fc  k), 

where  the  P,  is  that  described  in  Eq.  (A2-l). 

KJ 

Now  the  essence  of  the  pert\irbation  calciilation  is  to  assme  that  the  eigen¬ 
values  and  eigenvectors  for  the  perttirbed  matrix,  P  ,  in  Eq.  (A2-2),  differ  only 

„(0) 

slightly  from  those  of  P  J  and  that  they  can  be  represented  as  power  series 
expansions  in  the  small  parameter,  €,  with  zeroth-order  values  given  by  the  corres¬ 
ponding  forms  relative  to  I.e.,  if  the  m^^  eigenvector  and  eigenvalue  of  P 
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are  con5>onent,  and  respectively,  we  ass\me  that  we  can  write 


“ak" 


z  «=  z 
m  m 


^°K  ... .  f  .%<”) 

m  m  m  ^  ni 

n*0 


(A2-9) 


^  s  +  ...  =  Y.  €  + 

m  m  m  m  ^ 

n=0 


(For  the  usual  non-degenerate  case,  we  would  take  z^^^  *=  u  and  =7  :  we 

m  a  m  m 

aaipllfy  on  this  below.)  If  we  insert  Eq.  (A2-9)  into  Eq.  (A2-2),  recognize  a 
zeroth-order  identity  on  the  basis  of  Eq.  (A2-5a),  euid  collect  all  terms  in  like 
powers  of  €,  we  have  a  term-by-term  power  series  identity  relating  the  hi^er-order 
pertvirbation  correction  factors.  Confining  ourselves  to  the  first  two  orders,  we 
have 


First  order: 

'  m  '  m  m  m  m  m 


Second  order:  r<°>  *  P'  e  e  z<°> 

'  m  m  m  m  m  m  m  m 


(A2-10) 

(A2-11) 


At  this  point,  the  remainder  of  the  effort  consists  in  recognizing  and  using 
the  fact  that  z^^^  and  are  themselves  vectors  in  the  M-dimensional  space  spanned 

by  the  mutually  orthogonal,  normalized  set  >  and  hence  that  they  can  be  e:q)anded 

as  a  s\ai  over  the  with  appropriate  coefficients.  However,  first  we  must  con¬ 

sider  again  the  degeneracy  problem.  It  is  likely,  and  indeed  common  (e.g.,  the 
"line -splitting"  in  quantum-mechaniceJ.  problems)  that  under  perturbation  a  zeroth- 
order  degenerate  eigenvalue  will  give  rise  to  non-degenerate  (tonequal)  perturbed 
eigenvalues.  Now,  although  for  a  Q-fold  degeneracy  in  an  eigenvalue  of  a  Hermitian 
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matrix,  Q  mutually  orthogonal  eigenvectors  can  be  constructed,  any  linear  combina¬ 
tion  of  these  also  satisfies  the  eigenvalue  eqmtion  for  the  particular  degenerate 
eigenvalue,  and  may  substitute  as  an  eigenvector.  Thus,  although  we  know  that  the 
perturbed  eigenvector  has  as  its  zeroth-order  form  some  particular  such  linear  com¬ 
bination,  we  do  not  know  a  priori  what  this  form  is;  generally  it  will  not  corres¬ 
pond  to  a  particular  one  of  the  set  of  mutually  orthogonal  eigenvectors  which  we 
have  somewhat  arbitrarily  constructed  in  the  degenerate  zeroth-order  situation. 

Thus,  for  exanqple,  if  (say)  the  first  Q(m  =1,  ...  Q)  eigenvectors  of  describe 

such  a  degenerate  eigenstate,  we  must  replace  z^°^  in  Eqs.  (A2-9)  to  (A2-11),  for 
any  m  <  Q  (i«e.,  in  seeking  a  pert\arbed  state  corresponding  to  one  of  the  zeroth- 
order  Q-fold  degenerate  eigenstate)  by  some  linear  combination 


.(0)  , 


?  “  i  u. 

J=1 


(A2-12a) 


where  the  a^  will  also  have  to  be  determined  by  the  perturbation  calculation.  We 
still,  of  coiu-se,  take 


^m  =  ^m^=  ’'l  -  ’'2  =  •••  “  V 


(A2-i2b) 


Furthermore,  piursuing  our  earlier  remeurks,  we  expand  z^^^  and  z^^^  in  terms  of 
the  However,  since  we  have  alreeidy  introduced  in  the  zeroth-order  term,  those 

iij  which  correspond  to  the  •unperturbed  state,  we  argue^ ^ ^ ^ ^  that  we  can  omit  such 


We  ■will  focus  our  attention,  for  notational  convenience,  on  these  first  Q  eigenstates, 
where  Q  is  the  order  of  degeneracy  associated  ■with  the  first  eigen^value.  It  ■will  be 
tri-vially  ob^vious  to  the  reader  that  these  considerations  apply,  with  minor  notational 
changes,  to  calculating  perturbations  for  all  of  the  eigenstates. 
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terms  in  expanding  and  i.e.,  we  write 

m  m  ' 


,(1) 


m 


M 

I 


b  .  u. 

mj  i 


(A2-13) 


.(2) 


m 


M 

y'  c  u 

j4i 


(A2-14) 


(it  t\ims  out  that  if  one  attenrots  to  keep  terms  in  and  z^^^  relating  to  the 

mm 

zeroth-order  eigenstate,  the  associated  coefficients  are  basically  undetemined  by 
the  perturbation  calculation.  One  can  argue  that  their  value  is  then  si5)plied  by 
a  normalization  of  the  perturbed  eigenvector;  and  this  shows  them  to  be  always  higher- 
order  terms  which  never  enter  the  calculation  at  the  order  being  computed.  Although 
there  are  more  sophisticated  axqproaches  which  over-all  appear  to  deal  more  precisely 
with  this  problem,  for  the  low  order  to  which  the  pertxu'bation  corrections  would  v.r>ed 
in  our  case,  the  more  conventional  approach  being  presented  api>ears  quite  adequate 
and  purposeful.)  Further,  using  Eq.  (A2-5a),  it  now  follows  that 


t  b.r,u, 


m 


(A2-15) 


M 


r(o),u).  f 


'mj  ^5 


(A2-16) 


Thiis,  we  can  now  rewrite  Eqs.  (A2-10)  auid  (A2-ll)  as  follows; 


f-  b  ,  r,  u,  .  y  a  ,  r'  u,  .  £  b  ,  r  u,  +  ^  a  ,  u, 

«3  J  J  jti  U  J  mj  m  J  jtl  “J  ”  i 


(A2-17) 
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M  M  I  M  M  /  ^ 

y  C.7.U.  +  Y  b.i'u.  “  y  C.7U.  +  Y  b.i|f''u. 

J=Q+1  “J  J  J  j^l  “J  J  j-^1  oJ  “  J  oj  m  J 


+  y  ®-  A  (A2-18) 

A  oj  ’m  J  ' 


Finally,  ve  now  apply  the  orthononnality  of  the  set  >  i*e. 

It  is  also  convenient  to  note,  by  using  Eq.  (A2-6b),  that 


(A2-19) 


4  “j  -  P 


kj 


(A2-20) 


We  can  thus  obtain  from  Eqs.  (A2-17)  and.  (A2-l8)  (recalling  also  that  Y  ■  y  =  ...=  Y  )j 

X  2  ^ 


First  Ord.er: 


Second  Order: 


kiQ 

r  a  r*  «  a 
^mj  kj  mk^^m 

(A2-21) 

k>Q 

(A2-22} 

k<Q 

^  1  '  o) 

(A2-23) 

k>Q 

^mk*^  "  “mk^m  ^mk’^m 

(A2-24) 

Eq.  (A2-21)  may  be  regarded  as  a  set  of  Q  homogeneous  linear  algebraic  equations 
in  the  variable  a^j^jS 


t.  (  f'y  -  6,^  t^l))  -  0, 


J»1 


k  ■  1,  •  •  •  ,  Q 


(A2-25) 
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Non-trivlal  solutions  exist  only  for  values  of  such  that 

D 

(f’w  -  ♦i^’)  “  ■  ° 

where  the  subscript  Q  is  written  to  indicate  that  it  is  a  paxticular  submatrix  of  F 
which  is  involved  in  this  calculation.  The  Q  solutions  for  represent  the 

corresponding  first-order  correction  factors  for  the  initially  Q-fold  degenerate 
state.  If  the  values  for  are  unequal,  the  perturbed  state  is  non-degenerate 

if  some  are  still  equal,  these  conqsrise  a  residually  degenerate  state,  at  least  to 
first  order.  It  is  interesting  to  note  that  for  a  non-degenerate  zeroth-order  state 
(Q  *  1  in  Eq.  A2-26),  we  have  directly 

»  r  '  -  0  (A2-27) 

m  *  inni  ' 


(we  here  write  m  in  general,  rather  than  m  =  1  as  would  be  specifically  required  in 
Eq.  (A2-26)  by  our  predicate  that  we  are  examining  the  perturbations  of  the  first . 
ft-fold  degenerate,  state,  with  m  i.Q). 

For  a  given  solution  for  in  Eq.  (A2-26),  one  may  solve  Eqs.  (A2-25)  for 
the  corresponding  set  *  with  one  of  the  a^  left  undetermined.  If  there  is 

a  residual  Q*  -fold  degeneracy,  there  will  be  a  corresponding  set  of  q'  undeter¬ 
mined  coefficients;  the  Q*  euiditional  coefficients  can  be  determined  by  carrying  the 
computation  to  a  hi^er  order  (if  this  removes  the  degeneracy),  or  else  constructing 
an  arbitrary  mutually  orthogonal  set,  as  is  normally  done  in  connection  with  a 
degenerate  eigenstate.  In  any  case,  there  is  always  one  coefficient  left  unspeci¬ 
fied,  which  can  be  fixed  by  a  normalizing  requirement,  e.g.,  that  the  eigenvector 
z^^^  defined  in  Eq.  (A2-12)  be  normalized  to  unit  length. 


i 


1 


(A2-28) 
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We  can  next  determine  the  (h>Q)  from  Bq.  (A2-22), 


mk 


(A2-29) 


Thus  far,  ve  have  discussed  only  the  first  order  perturbation.  We  can  now 
return  to  the  second  order.  From  Eq.  (A2-23),  if  the  m^^  state  is  one  for  which 
the  degeneracy  was  removed  in  the  first  order  (hence  the  a^  determined) 

(any  k  <  Q  may  be  used  here).  If  the  m^^  state  is  still,  in  first-order,  one  of  a 

o' -fold  l^rrenerate  system,  we  have  to  InstoT-rl  eii^loy  the  solution  for  the  b  in  torina 

n, 

t 

of  the  Q  \andetermlned  as  given  by  Eq.  (A2-29);  we  substitute  this  into  the 

appropriate  q'  equations  of  Eq.  (A2-23)  and  end  up  basically  with  a  set  of  Q* 
homogeneous  linear  algebraic  equations  in  Q*  variables.  The  requirement  that 

the  determinant  of  this  Q'-fold  system  vanish,  for  non-trivial  solutions  to  exist, 

(2) 

then  gives  the  values  for  ,  and  then  in  turn  the  equations  yield  the  correspond¬ 
ing  sets  of  Sjuj,  with  one  vindetermined  and  to  be  fixed  by  normalization.  If  there 
is  residual  degeneracy  in  this  order,  it  must  be  handled  Just  as  commented  vq>on 
earlier  (either  look  for  higher  order  terms,  or  settle  for  the  degeneracy  and  con¬ 
struct  suitable  a^  accordingly). 

The  determination  of  the  c^  from  Eq.  (A2-2U)  is  then  straightforward.  By  way 
of  further  concluding  our  resvilts,  we  note  that  wherever  in  our  perturbation  formulas 
a  factor  appears,  a  factor  C  appears  in  conjunction  with  it.  Thus,  we  see  that 

jK 

exactly  the  forms  in  Eq.  (A2-8b)  always  enter  in.  Thus,  for  instance,  we  note  that 
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in  the  first  order  pertiirbation, 


m 


z('». 

m  n 


Q 

I 

J=1 


a  .  u.  + 

mj  J 


M 

I 

i^i 


Q 

^  JP 


u 


J 


Or 


m 


a  .  u.  + 
i 


M 

z 


p?i  %  ^jp  s 


(A2-30) 


Similar  remarks  hold  for  the  eigenvalue,  and  into  all  higher-order  pertvirbations . 

Finally,  ve  again  ceLll  to  attention  that  ve  have  basically  described  the  con- 
stinction  of  one  perturbed  eigenstate,  or  of  one  member  of  a  perturbed  degenerate 
eigenstate.  The  calculation  must  be  carried  in  turn  out  for  each  of  the  total  of  M 
eigenstates  associated  with  the  f -matrix. 
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APPENDIX  3:  OPTIMIZATIQN  OF  A  7ZPICAL  ANGLE  DIVERSITY  ARRAY 
a)  Radiation  Efficiencies 

We  will  consider  as  a  typical  angle  diversity  receiving  array  that  pictured  in 
Figure  5«  There  is  a  central  heam  (Beeun  3)  and  four  symmetrically  placed  side 
beams  (Beams  1,  2,  4,  5)*  On  grounds  of  physical  symmetry,  we  will  assume  the 
latter  four  are  all  a  priori  assigned  equal  radiation  efficiencies, 

(A3-la) 

while  the  center  beam  has  some  different  efficiency,  by  a  ratio,  K, 

=  Kq  (A3-lb) 

Furthermore,  we  assxime  all  the  beams  to  be  excited  by  currents  in  a  common  aper- 
txire,  with  only  phasing  across  the  aperture  to  produce  the  different  beam  point- 
Ing  directions.  In  this  case,  we  can  take  all  the  beam  patterns,  4o  be 

real- valued,  and  all  Identical  in  form. except  for  angle  of  pointing.  We  also  assume 
that  al2  have  one  common  lineai  polarization.  In  this  case,  bj  the  symmetry,  we  can 
write 

Pll  “  ^22  “  ^33  “  ^44  “  ^55  “  ^ 

Pl2  “  P25  “  ^5  ‘  ^14  “  “1 

(A3-2) 

^13  “  ^23  “  ^3  “  ^53  “  “2 
^5  "  ^24  “  ®3 

For  the  numerical  example,  we  will  consider  a  case  where 

1  >  Og  >  Cl  >  03  >  0  (A3-3) 
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On  physical  grounds,  such  an  inequality  will  usually  hold  for  the  absolute  values 
of  the  Ofj.  However,  as  noted  earlier  in  Appendix  1  (Figure  3)^  some  of  the  0£j  may 
be  negative,  and  indeed,  the  inequality  may  not  8Q.vay5  hold,  even  with  the  configu¬ 
ration  assumed.  However,  it  would  more  usually  be  the  case  than  not  with  this  con¬ 
figuration,  aind  fortunately,  the  validity  of  Eq.  (A3-3)  for  o\ir  nianerical  values 
serves  well  to  simplify  the  later  arithmetic  in  our  Illustrative  example. 

Then,  along  with  Eqs.  {A3-1),  we  have 


The  eigenvalues,  of  P  are  determined  by  solving  det(  P  -Yl)  “  0 


(A3-5) 
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This  may  be  slngjllfied  by  factoring  K  {K  ^  o)  from  the  third  column,  and  K 
from  the  third  rowj  and  also  factoring  |q|^  from  every  row.  Defining 


we  then  have 


1-x 

“1 

“2 

“1 

“3 

“1 

1-x 

“2 

“3 

“1 

“2 

“2 

X 

“2 

“2 

A.  “ 

CM 

“1 

“3 

“2 

1-x 

“1 

a. 

a. 

Oo 

a. 

1-x 

(A3-6) 


(A3-7) 


We  note  at  this  point  that  K  only  enters  throu^  |k|^  and  hence  with  no  loss  in 
generality,  we  may  hereafter  regard  K  as  real.  Next,  Eq.  (A3-7)  can  be  simplified 
by  subtracting  the  last  column  from  the  first,  and  the  fourth  from  the  second: 


0  *» 


1-a^-x 

0 

«2 

“1 

“3 

0 

l-a^-x 

«2 

“3 

“1 

0 

0 

1  - 

“2 

“2 

0 

a^+x-1 

«2 

1-x 

“1 

ct^+x-l 

0 

“2 

“1 

1-x 

(A3-8) 
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Next,  adding  the  first  row  to  the  last,  and  the  second  to  the  fourth: 


l-a^-x 


l-a^-x 


a_ 


1  -  -2^ 


2a„ 


a. 


a. 


20^ 


a. 


a. 


a_ 


14a2“X  20^ 


20£^  l-*ci^-x 


Hence,  ve  can  Immediately  reduce  the  equation  to 


0  -  (l-a.-x)‘ 


1  -  -2S- 


ct« 


20t2  I4a^-x 


2Q!„ 


ot. 


2a, 


20^  l+a^-x 


(A3-9) 


If  we  now  subtract  the  third,  row  from  the  second,  and  then  add  the  new  second 
column  to  the  third. 


0  ■  (l-a^-x)' 


1-4 


2a„ 


2a„ 


l*Cl^-2PC^~x 


^1  l+2a^+a2-x 


(l-a2-x)^(  I4a2-2a^-x) 


1  -  ^ 
r\ 


2a„ 


20^ 


l+20^4a2-x 
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Thus  the  five  roots  sore: 


■  Xg  =  l-O!^  (double  root) 


x^  « 


and  x^,  x^  are  the  solutions  to 


(1  -  (1+20^40!  -x)  -  l40^  -  0 


x^-x(K^+l+2a^+a  )  +  Y?  (i+2a^+a  )  -  4l^  -  0 


(A3-12a) 


‘4^  (K^-n-f2a^4a^)  jh  y[^-m-2a^+a^' ^  4  4k^  [40^  -  (n-2a^+a^)] 


vhlch  can  alternatively  be  vritten 


^4 1  ^  Vf  -  (1+20^40^)]^  4  16K^  0^ 


(A3-12b) 


We  note  that  in  Eq.  (A3-I2)  the  leurger  eigenvalue  is  Xj^  and  the  snsiller  is  x^. 
From  the  latter,  since  on  physical  grounds  it  cannot  be  negative,  ve  see  that  ve 
mvist  have 


4o^  -  (l+2a^4a  )  <  0 


p  1+20  +a 
CC  <  - r= — ^ 

2  4 


(A3-13) 


Althou^  not  obvlo\is  by  Inspection,  this  latter  inequality  probably  follows  fron 
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the  haelc  aBBvut^tlons  of  our  configuration. 

Finally,  now,  from  Eq.  (A3-12b),  we  note  that 

Or,  for 

<  1+20^40^  ,  xj^  >  l+aa^+a^ 

>  1+20^40^  ,  x,^  >  >  l4-2a^40£2 

In  either  case  then, 

Xj^  >  14-20^402  >  x^  or  x^ 


(K^4-i4-2a^4a  )  4-  |k^  -  (l4-2a^4a  ) 


(A3-14) 


and  hence  Xj^  is  always  the  largest  eigenvalue. 


Recalling  Eq.  (A3-6)  and  Eq.  (A3-9)>  ve  thus  have  Insediately 


(>'kW 
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(It  is  incidentally  readily  verified  that  |q|^  <  1  also,  as  is  required 
physically. ) 

Let  us  now  assume  the  equality  sign  in  Eqs.  (A3-I5)  and  (A3-16),  i.e.,  that 
whatever  the  value  of  |k|,  we  have  a  lossless  antenna  system  so  designed  that  there 
are  no  cross-covqplings  or  losses  other  than  fundamentally  inherent  and  implied  by 
the  beam  overlaps  and  the  conservation  of  energy  limitations . 


b)  Optimization  of  Diversity  Operation 

We  now  refer  to  the  optimization  problem  discussed  in  Section  8.  We  will 
assiime  that  the  matrix,  G,  of  beam  correlation  factors,  G.  (defined  in  Eq.  50)  has 
the  same  kind  of  symmetries  that  the  S-matrlx  had.  This  is  basically  the  assim^- 
tlon  that  g(Q)  in  Eq.  (50)  is  circularly  symmetric,  and  that  Beam  3  is  pointed  at 
its  center  of  symmetry.  This  assvm5)tion  is  not  a  particularly  accurate  descrip¬ 
tion  of  the  tropospheric  scatter  fields,  which  display  a  different  scale  of  varia¬ 
tions  in  the  veirtical  than  in  the  horizontal  directions.  However,  since  our  pur¬ 
poses  here  are  purely  illustrative,  we  believe  the  slnpllflcatlon  in  resulting 
mathematics  Justifies  the  symmetry  assuaiptlon  on  g(0).  The  point  is  that  with 

this  assumed  symmetry,  we  can  immediately  note  symmetries  in  G  .  which  are  almost, 

JK 

but  not  coopletely  anaJ.ogous  to  those  in  Eq.  (A3-2)}  we  write  these  as 

If 

=  h 
33  o 

®11  =  °22  ^  ^ 


°12  =  °25 

^13  “  °23  “  ®43  ■'^53*^ 


(A3-17) 


G15  =  Ggi^  =  h3 


A3-7 


ARM-237 


Hence,  we  can  determine  the  eigenvalues,  for  the  L-matrix  (Eq.  51)  in 


eOinost  direct  analogy  to  those  determined  for  P  •  Thus,  if  we  define 


o 


(A3-l8a) 


*"0 


^o  1^ 


(A3-l8b) 


we  obtain  as  the  solutions  for  the  eigenvalues 


A '  4  ■  "o  |i|^  ■  ”0  i'jp  ^2  ■  “o 


(A3-19a) 


(A3-19b) 


2  (  K’  ^+l+2h^+h- )  +  y  (  k'  ^+l+2h^+h-  )^+4k'  ^  Uh|-(H-2h^th^) 


(A3-19c) 
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p 

In  these,  jq|  is  now  assumed  to  be  given  by  the  equality  sign  in  Eq.  (A3-I5) 
(in  terms  of  K,  not  yet  specified).  Furthermore,  the  error  rate  function  F(z) 
to  be  minimized  for  optimum  performance  (Eq.  47),  is 


“  2  detU+zL) 


(A3-20) 


(z  will  be  further  specified  below.)  With  L  Hermitian,  one  knows  that  a  diagonaliz¬ 
ing  unitary  matrix  exists  for  L,  say  T,  such  that 


=  I 

and 

T'^LT  «  Ji 

where  JLi.s  the  diagonal  matrix 


But  then,  since  det  T  =  1, 

det(l+zL)  =  det(T'^IT  +  zT'^ET)  •=  det(l+z/) 


(A3-21) 


(A3-22) 


Hence, 


(A3-23) 


(A3-24) 


=  noise  level  in  individual  receiver  (same  for  all  receivers) 
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and  ve  have  F(z)  as  an  e:q>llclt  function  of  K  throvigh  Eqs.  (A3-I9)  and  (A3-13). 

To  naximlze  the  res\iltlng  e^qpresslon  algebraically  Is  tedious.  Therefore,  ve  nay 

take  vhat  may  be  regarded  as  a  typical  set  of  nmbers  for  (see  Appen- 

"  *  »  »  « 

dlx  1),  and  another  typical  set  for  h^,  h^,  hg,  h”.  With  these,  wt  wil]  plot 
in  Flguree  7  and  3  the  variation  of  h'(z'  xs  s  function  of  K  foi  variouE  assumed 
average  signal-noise  ratios. 

The  particular  numerical  values  taken  are  for  gaussiem  beam  patterns,  in  the 
configuration  of  Figure  3)  where  Beam  3  intersects  the  other  beams  at  their  33dlb 
points.  Prom  Figures  2  and  3,  ve  then  obtain 

=  O.U9 

Og  «  0.2U 

=  0.06 

For  the  intensity  pattern  g{0),  we  assume  a  form 


where  0,  normalized  to  the  same  angle  coordinates  used  in  Figures  2  and  3;  ^ 

effect  the  source  "spread"  normalized  to  the  beam  width  of  our  individual  beams. 

The  value  determines  the  actual  power  level  of  the  received  signals.  The  values 
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obtained  for  the  G  matrix  are  then: 


e  -  1 

o 

0-2 

0 

h"  =  0.1574 
o  o 

h"  =  0.2697  Q 

h'  =  0.0644  Q 

0  0 

h'  -  0.1832  Q 

0  0 

«  0.0442 

h^  =  0.0984 

h^  =  0.0833 

h^  =  0.1632 

h^  -  0.0306 

h^  -  0.0524 

The  Individual  vUl  thus  be  proportional  to  Q^.  Thus,  for  each  term  In 
Eq.  (A3-2U), 

zl  ."ts-i' 

*■^1  2N  1 
o 

where  J' ^  is  the  value  of  the  eigenvalues  if  »  1  is  taken  in  defining  the  elements 

"  ^o 

of  the  G-matrix.  Clearly,  furthermore,  h^  ^  is  the  signal-noise  ratio  which  would 
appear  in  the  receiver  associated  with  the  central  beam  if  ^  other  beams  were  pre¬ 
sent  (and  |k^|  |q|^  =  1  were  involved).  Thus,  let  us  redefine  s  as  the  latter  signal- 
noise  ratio,  and  rewrite  Eq.  (A3-25)  in  the  form 


(A3-25a) 


(A3-25b) 
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In  Figures  6-8  have  been  plotted  the  radiation  efficiencies  and  F(z)  as 
I  12 

functions  of  K  .  The  latter  variable  has  been  graphed  logarithmically,  to  bring 

out  the  essential  features.  In  almost  all  cases,  the  values  plotted  for  |Kp  >  0.1 

are  relatively  close  to  those  for  jKp  =•  0  (the  latter  have  been  indicatcu  via 

arrows).  In  Figure  6  are  given  the  reidiation  efficiencies,  |q|^,  for  the  four  outer 

beams,  and  |K|^|qj^  for  the  central  beam.  It  is  clear  that  for  |k|^ - »0,  the 

central  beam  is  effectively  omitted,  and  the  residual  limitation  on  |q|^  is  the 

I  12 

Interaction  only  among  the  four  outer  beams.  Likewise,  for  K  very  large, 

|k|^  jq|^ - >1  while  |q|^ - >0,  and  in  effect  only  the  central  beam  is  present.  In 

Figures  7  and  8  are  given  plots  of  the  error  rate,  ^(z),  for  selected  values  of 
z  (O,  3,  10,  15>  and  20  db),  and  for  6^  =  1  and  0^  =;  2,  respectively.  The  ordinate 
is  broken,  since  for  the  range  of  z  taken,  the  error  rates  vary  over  several  decades, 
whereas  it  seems  desirable  to  plot  all  the  curves  on  one  graph. 

The  most  striking  observations  from  Figures  7  and  8  are  that  there  is  an  opti- 

1  i2  1  |2 

mum  value  for  |K|  ,  generally  around  |  K|  =  0.7>  and  somewhat  independent  of  z  or 
0^^  that  the  error  rate  ctirves  have  a  broad  minimum  around  the  optimum,  so  that  the 
choice  is  not  criticalj  that,  indeed,  for  the  configuration  assumed,  and  the  signal- 
noise  ratios  examined,  the  variation  of  error  rate  is  by  no  more  than  a  factor  of  2 
over  the  range  0.1  <  5  Changes  in  error  rate  by  a  factor  of  2  are  equivalent, 

due  to  the  exponential  dependence  of  error  rate  on  signal-noise  ratio,  to  very  minor 
changes  in  signal-noise  ratio.  Thus  one  may  especially  notice  that  in  the  configvura- 
tion  examined,  the  optirman  presence  of  the  center  beam  (optimum  |k|  )  is  very  little 
better  than  "Che  complete  absence  cf  the  center  beam  |k|^  ==  O).  That  is,  the  added 
losses  (due  to  cross-coupling)  in  the  power  received  with  the  outer  beams,  when  the 
center  beam  is  active,  alra'''St  conpletely  coitpensate  for  the  additional  diversity 
advantage  of  adding  this  center  beam,  (indeed  for  the  0^  =  2  case,  the  optimum 
IkI""  for  small  signal-noise  ratio  is  in  fact  |k|^  =  0.) 
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It  is  known  that  for  sufficiently  low  error  rates  (sufficiently  high  signaJ.- 
noise  ratio),  the  added  diversity  advantage  of  having  five  rather  than  four  beams 
must  appear.  This  is  indicated  by  the  decreasing  width  of  the  minimum  as  signal- 
noise  ratio  Increases.  However,  for  our  exaaiple,  the  error-rate  at  which  this 
diversity  advantage  appears  is  apparently  far  below  those  which  appear  in  our 
curves,  which  in  turn  are  typical  of  operational  requirements. 

One  may  also  question  whether  the  beam  overlap  factors  used  here  are  realistic. 
We  have  not,  in  all  this,  specifically  derived  the  related  cross -covqpllngs  between 
the  feed-lines.  In  an  actual  design,  as  mentioned  earlier,  the  design  engineer 
might  attenq)t  to  modify  the  feed  structure  to  reduce  such  cross-co\qollngs,  in  a 
way  which  woiiLd  not  alter  the  main  lobe  pattern  but  might  drasticedly  alter  the 
side-lobes  and  hence  the  overlap  factors.  (One  may  also  question,  on  the  other 
hand,  whether  this  is  really  ever  done  in  designing  an  eurray  for  receiving  only; 
the  penalties  are  obvious  for  a  transmitting  array,  but  it  is  not  apparent  that 
they  have  ever  been  regarded  seriously  for  receiving-only  arrays . ) 

At  any  rate,  this  exaiiq>le  sho\ild  be  taken  only  as  an  illustration  of  the  kind 
of  effects  which  have  been  ln^uted  from  the  mathematics,  and  as  an  indication  of 
the  kind  of  detailed  investigation  idiich  should  be  undertaken  in  connection  with 
any  angle-diversity  designs. 
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ON  CROSS-COUPLING  IN  MULTIPLE  BEAM  ANTENNAS 

PART  II  CHARACTERISTICS  OF  THE  CROSS-COUPLING  MATRICES 


by 


S.  Stain  and  J.E.  Storer 


ARM238 


ABSTRACT 

The  analysis  in  Part  I  has  been  extended  to  the  characterization 
of  the  scattering  matrix  for  a  multiple -beam  antenna  system.  The 
constraints  implied  by  beam  overlap,  and  the  principle  of  conservation 
of  energy,  are  derived.  The  more  specific  properties  are  given  for 
antennas  further  characterized  as  lossless,  reciprocal  systems.  Some 
directions  for  further  useful  research  are  pointed  out. 
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1 .  Summary  of  Prior  Res\ilts 

In  Part  I  ,  the  principle  of  conservation  of  energy  was  applied  to  a  multiple 
beam  antenna  (M  beams).  An  Hermitian  positive  semi-definite  MxM  matrix,  p,  of 

beam-overlap  factors  was  defined,  and  a  set  of  M  beam  excitation  factors,  , 

2  th  <■ 

for  which  I  I  represent  the  radiation  efficiency  for  the  k  beam.  Another  MxM 

Hermitian  positive  semi -definite  matrix,  F,  was  defined,  with  its  elements  given 

by 


^kj  \  *^kj 


(1) 


Furthermore,  the  antenna  system,  considered  as  a  junction  among  the  M  feed-lines 
(each  excites  exactly  one  of  the  beams),  is  characterized  by  an  MxM  "scattering- 
matrix’^  S,  whose  elements  S  .  represent  the  mutual  coupling  coefficients  (and 
self-reflection  coefficients)  among  the  set  of  feedlines.  For  a  given  set  of  incident 
waves  in  the  several  feed-lines,  with  amplitudes  ^x^J-  characterized  as  a  column 
vector,  X,  the  set  of  reflected  waves  is  denoted  by  another  set^y^^,  or  an 
equivalent  column  vector,  y. 


y-Sx  (2) 

The  speciiic  result  obtained  from  applying  the  principle  of  conservation  of  energy 
involves  the  eigenvalues  and  eigenvectors  of  F.  Since  F  is  Hermitian,  and  positive 
semi-definite,  its  M  eigenvalues  are  known  to  be  positive  real  or  zero;  we  denote 
by  y  the  diagonal  matrix  (i.  e. ,  matrix  with  non-vanishing  elements  only  on  the 
main  diagonal)  whose  diagonal  elements  are  these  eigenvalues  of  F, 


S.  Stein,  ''On  Cross-Coupling  in  Multiple-Beam  Antennas,  Parti,  Radiation 
Efficiencies  (With  Applications  to  Angle  Diversity),  "  ARM  No.  237,  3  March  1961. 
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The  y^,  >  0  the  aforementioned  eigenvalues.  Furthermore,  the  eigenvectors 

of  r  form  a  mutually  orthogonal  set  (or,  for  a  degenerate  eigenvalue,  can  be 
constructed  as  such),  and  can  further  be  normalized  to  unit  length.  Within  a  lack 
of  uniqueness  discussed  in  Appendix  1,  one  can  define  a  unitary  matrix,  U,  whose 
columns  are  the  respective  orthonormalized  eigenvectors  of  F,  and  which  has  the 
properties 


U'^U  =  I 
u'^ru  =y 
r  =  UyU"^ 


(I  is  the  identity  matrix,  I..  =6..) 


(4) 


The  fundamental  equation  derived  in  Part  I  (Part  I,  Eq.  36)  was  the  inequality 


z'^(I  -  y)z  >  z'*’(U’^S’^SU)z  (5) 

where  the  equality  sign  holds  if  and  only  if  a  lossless  antenna  system  is  being 
considered.  The  immediate  deduction  from  Eq.  5  was  that  all  the  eigenvalues  of  F 
are  restricted  by  the  inequalities 


Most  of  Part  I  was  devoted  to  exploring  the  consequences  of  Eq.  6,  with  respect  to 
limitations  imposed  by  beam-overlaps  (by  the  form  of  (3)  upon,  the  achievable 
radiation  efficiencies.  Once  the  latter  are  chosen,  consistent  with  Eq.  6,  the 
are  completely  specified  within  the  assumed  beam  patterns. 


Again  we  use  the 

:!< 

A 


notations,  for  any  complex  matrix  A, 

❖  sit 

=  complex  conjugate  of  A;  (A  ). .  =  A  .. 

^  J 


=  transpose  of  A;  (A  )^j  =  A^^ 

=  conjugate  transpose  of  A;  (A  )^j 
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In  this  report,  we  will  explore  the  further  implications  of  Eqs.  5  and  6, 
with  respect  to  characterizations  of  the  scattering  matrix.  We  will  provide  a 
rather  complete  characterization  for  the  very  important,  most  common  case  of  a 
lossless  reciprocal  antenna  system,  which  includes  all  antenna  systems  comprised 
of  essentially  passive  lossless  elements.  For  more  general  systems,  we  will  only 
be  able  to  indicate  a  partial  characterization,  which  includes  all  the  implications  of 
conserva.tion  of  energy  and  the  beam  specifications,  but  apparently  always  requires 
other  constraints  (such  as  losslessness,  reciprocity)  to  complete  the  description 
of  the  system. 


2.  General  Solution  for  S 

In  Eq,  5  (along  with  Eq.  6),  we  recognize  the  I-y  is  a  diagonal  matrix  with 
elements  positive  real  or  zero.  We  can  thus  define  a  diagonal  matrix,  A,  whose 
elements  are  the  set  of  positive,  real  or  zero  numbers 


k  =  1,  .  .  . ,  M 


(7) 


A,  ,  6,  , 


kj 


kj 


0  \ 


A  - 


•  X 


M 


/ 


We  note 


T  ’i'  + 

A  A  A  ^  A^ 


v7e  can  then  write 
I  -  y  -•  A^ 


and  Eq.  5  becomes 


z'^A^z  >  z'^(U'^S‘*'sU)z 


(8a) 


(8b) 


(8c) 


(9) 


(10) 


Now  in  Part  I,  it  was  shown  that  with  all  other  relevant  parameters  fixed, 
the  maximum  radiation  efficiencies  for  all  the  beams  are  attained  when  the  largest 
of  the  is  equal  to  unity  (the  largest  value  com.patible  with  Eq.  6).  Both  the 
diagonal  matrices  (I-y)  and  A  would  then  have  some  of  their  diagonal  elements 
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zero  (as  many  such  elements  as  the  order  of  degeneracy  of  the  largest  eigenvalue). 
To  avoid  the  mathematical  complications  necessitated  by  accounting  for  this  effect, 
we  will  assume  in  this  section  that  none  of  the  X,  vanish  (all  y.  are  less  than  unity), 

r. 

On  this  basis,  A  has  a  well-defined  inverse  A  ,  which  we  will  use.  However  this 
inverse  will  be  seen  to  be  absent  from  the  final  result,  and  one  will  be  able  to 
argue  then  that  the  result  applies  as  well  to  the  limiting  case  where  some  of  the 
X^  do  vanish.  In  Appendix  2,  we  further  indicate  by  paralleling  the  derivation  in 
this  section,  a  more  rigorous  demonstration  of  this  assertion. 

Now,  assuming  A  ^  exists,  let  us  define  a  new  matrix 


T  =  SU  A  (11a) 

S  =  TAU'^  (11b) 

so  that  Eq.  10  becomes 

z'^A^z^  z'^A  T'^TAz  (12) 

Defining  a  new  (equally  arbitrary)  vector 

£  -  Az  (13) 

we  can  further  write  the  inequality  as 

£^£>£^T^T£  (14) 

Clearly,  the  matrix  T  T  is  Hermitian,  and  positive  definite,  since  S^S  is  positive 
definite  (shown  in  Part  I).  Thus  it  can  always  be  written  in  the  form 

t'^T  =  (15a) 

where  P  is  an  appropriate  unitary  matrix 

P'^P  =  1  (15b) 


and  p  is  a  diagonal  matrix  whose  diagonal  elements,  are  positive  real  or 

zero.  The  requirement  imposed  by  Eq.  14  is  then  readily  observed  (see,  e.  g,  , 
Appendix  2)  to  be  simply  that 

^  1  for  all  k  (16) 

Since  the  latter  statement  is  the  totality  of  information  contained  in  the  conservation 
of  energy  equation,  Eq  14,  we  see  that  the  class  of  all  permissible  solutions  for 
T  T  is  hencq  given  by  Eq.  15,  where  P  is  any  unitary  matrix,  and  any  diagonal 
matrix  satisfying  Eq.  1  6, 
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It  immediately  follows  by  ’’factoring"  Eq.  15a  that  within  this  class  of 
solutions,  T  itself  has  as  its  most  general  form 

T  ::  WmP  (17) 

where  W  is  again  any  unitary  matrix 

W'^W  =  I  (18) 

Using  Eq,  11,  the  corresponding  general  class  of  solutions  for  S,  consistent  with 
conservation  of  energy,  is  just 

S  =  WmPAU'^  (19) 

In  the  general  class  of  solutions  just  given,  only  the  diagonal  matrix  A  is 
completely  specified  (by  our  choice  of  the  multiple  beams,  and  the  choice  of 
radiation  efficiencies),  and  the  matrix  U  partially  specified  (see  Appendix  1). 

We  can  expect,  and  will  discuss  below,  that  the  arbitrariness  of  choice  of  the 
unitary  matrices  P  and  W  has  physical  implications  concerning  certain  freedoms 
of  choice  in  designing  the  antenna  system  to  drive  the  specified  set  of  beams.  The 
inequality  which  specifies  )x,  however,  is  a  source  of  difficulty.  In  the  lossless 
case,  of  course,  Eq,  16  is  to  be  taken  with  the  equality  sign,  and  we  will  see 
below  that  this  greatly  reduces  the  generality  of  possible  forms  for  S. 

In  the  more  general  lossy  case,  the  inequality  signs  in  our  equations,  Eq.  5 
for  example,  could  similarly  be  replaced  by  equality  signs  if  a  term  were  added 
which  accounted  for  all  losses  within  the  antenna  system  (considered  as  a  junction), 
i,  e,  ,  Eq.  5  rewritten  as 

z^(I  -  7)2  -  z^(U^S^SU)z  +  L  (20) 

where  L  is  a  term  accounting  for  all  the  losses.  If,  in  fact,  such  losses  were  all 
ohmic,  E  would  be  a  Hermitian  quadratic  form  in  z, 

L  ^  z'^(U'^L  U)z  (21) 

o 

(where  is  the  Hermitian  ma<^rix  of  the  quadratic  form  giving  the  losses,  when 
the  variable  is  the  set  of  incident  waves  in  the  feedlines).  If  L,  or  L^,  were  known, 
one  could  presumably  make  corresponding  specific  statements  about  p,  much  like 
the  ones  to  be  made  below  in  the  lossless  case.  If  the  losses  were  due  to  specific 
"pads''  inserted  into  the  networl^,  one  might  be  readily  able  to  describe  such  a 
matrix,  under  more  general  lossy  conditions,  however,  it  is  not  apparent  that 
enough  is  generallv  known  about  the  quadratic  form  to  make  use  of  it.  Results  in 
the  remainder  of  this  paper  will  be  confined  to  the  lossless  case. 
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Other  specializations  can  be  made,  not  related  to  the  lossless  vs.  lossy 
issue,  for  example  the  assumption  of  reciprocity.  Combined  with  losslessness, 
below,  this  yields  a  particularly  interestingdelineation  on  the  form  of  S.  The 
characteriz’.ation  assuming  reciprocity  only,  however,  does  not  appear  to  be  useful 
and  again  v/ill  not  be  cited  explicitly.  As  we  comment  in  our  conclusions,  there 
thus  seem  to  be  a  number  of  interesting  problems  yet  to  be  explored  in  this  general 
direction  of  system  constraints.  These  we  necessarily  regard  as  beyond  the  scope 
of  our  currently  achievable  results  presented  in  this  paper. 

4.  Lossless  Reciprocal  Systems 

By  far,  the  widest  class  of  antenna  systems  can  be  regarded  as  possessing 
both  the  properties  of  reciprocity  and  losslessness. 

For  a  lossless  system,  the  equality  sign  holds  in  Eq.  16,  and  we  have 

=  I  (22) 

Since  p  p^,  p  itself  can  be  replaced  by  any  unitary  matrix;  the  resulting  product 
WpP  in  Eq,  17  becomes  simply  the  product  of  three  arbitrary  unitary  matrices,  and 
leads  to  the  simple  statement 

T  -  any  unitary  matrix  (23) 

(This  result  is  also  apparent  by  simply  setting  p^  =  I  in  Eq.  15),  Thus,  the  property 
of  losslessness  reduces  to  a  statement 

S  =  QAU'*'  (24a) 

where  Q  is  any  unitary  matrix, 

Q^Q  -  I  (24b) 

In  order  for  the  network  to  be  also  reciprocal,  we  require 
S  =  (25) 

(note,  since  S  is  complex,  this  is  not  Hermitian  symmetry) 

Clearly  arbitrary  choices  of  Q  will  not  satisfy  this  requirement,  and  hence  would 
generally  imply  non-reciprocal  networks.  (Aside  from  the  difficulties  in  specifying 
particular  kinds  of  non-reciprocity,  if  one  may  use  such  vague  terminology,  it  was 
also  pointed  out  in  Part  I  that  such  systems  are  basically  irrelevant  in  the  important 
problem  of  analyzing  receiving  systems  in  terms  of  equivalent  transmitting 
properties,  )  Now,  to  satisfy  the  reciprocity  requirement  in  Eq.  25,  the  class  of 
solutions  in  Eq,  24  must  also  obviously  satisfy  the  relation  ^ 
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+  +  T  *  T 

QAU  =  (QAU  )  =  U  AQ^ 


(26) 


If  we  define  a  new  matrix,  F,  by 
Q  =  U''f 

then  F  must  also  be  unitary, 

F'^F  =  I 

and  satisfy,  from  Eq.  26, 

FA  =  Af”^ 

In  ierms  of  components,  since  A  is  diagonal,  Eq.  28b  becomes 

F.  .(X.  -  X.)  =  0 
ij  1  J 


(27) 


(28a) 


(28b) 


(29) 


I.  e.  ,  either  X.  =  X.  or  F. .  =  0.  Thus  if  all  the  diagonal  elements  are  unequal,  F 
^  J  ^  J 

will  itself  be  any  diagonal  unitary  matrix.  In  the  more  general  case,  if  some  of 
the  X^  have  equal  values  (corresponding  to  degenerate  eigenvalues  of  F),  off- 
diagonal  elements  are  possible,  at  intersections  of  the  corresponding  rows  and 
columns.  If,  without  loss  in  generality,  we  assume  the  X^  for  the  moment  to  be 
ordered  (grouped),  then  clearly  the  non-vanishing  elements  allowed  by  Eq.  29  must 
be  clustered  such  that  F  will  be  a  partitioned  diagonal  form. 


F  = 


(30) 


\ 


0 


where  the  size  of  the  aratiojiu,  F^,  is  equal  to  the  order  of  multiplicity  of  the 
corresponding  eigenvalues.  Within  this  context,  Eqs.  28  simply  now  require  that 
the  individual  partitions,  F^,  be  unitary  and  syn^metric  matrices  (within  their 
dimensionality) . 

However,  it  is  well  known  that  any  cy..,,.  ctric  matrix  can  always  be  factored, 
in  an  infinite  number  of  ways,  as  the  product  of  some  matrix  am.'  its  transpose, 
i.  e.  ,  for  each  F^  we  can  write 


F.  =  V.V. 
1  11 


(31) 
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Furthermore,  to  satisfy  the  unitary  property  of  F^,  it  is  sufficient  to  choose  only 
those  factorizations,  V^,  which  are  themselves  unitary 

V^V.  =  1  (32) 

1  1 

(it  is  readily  shown  that  such  exist).  Making  any  such  factorization  for  the 
F^,  we  can  then  obviously  write 

F  ^  VV'^  (33) 

where  V  is  a  diagonal  partitioned  unitary  matrix  of  form  similar  to  F, 


Then  S  has  the  form 

^  «L  ..L 

S  =  QAU  =  U  FAU  =  U  VV  AU  (35) 

But,  finally,  any  diagonal  partitioned  matrix  such  as  V  commutes  with  A  (analogous 
to  Eq.  Al-8)  and  thus  we  can  always  write  this  S  in  the  form 

S  =  uVaV^u"^  (36) 

=!«  >!<  T 

=  (U  V)A(U  V) 

But  (Appendix  1),  with  V  unitary  and  partitioned  as  shown,  U  V  is  simply  where 
is  some  other  unitary  matrix  which  diagonalizes  F;  i.e.  ,  within  the  lack  of 
uniqueness  of  the  matrix  U,  our  most  general  solution  for  S  in  the  lossless  re¬ 
ciprocal  case  has  the  form 

S  =  u’^AU'*'  (37) 


5.  Physical  Significance  of  the  Lossless,  Reciprocal  Characterization 

For  a  given  multiple -beam  design,  and  decision  on  beam  radiation  efficiencies, 
the  matrix  A  is  completely  fixed,  and  the  unitary  matrices,  U,  which  enter  Eq,  37 
are  fixed  to  within  certain  qua  si -diagonal  forms  (Eq.  Al-7),  It  is  useful  to 
interpret  the  meaning  of  this  arbitrariness  still  inherent  in  the  specification  of 
S  via  Eq,  37, 
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We  recall,  Eq.  2,  that  S  relates  the  set  of  reflected  waves  in  the  feed-lines 
to  the  set  of  incident  waves, 

y  =  Sx 

Using  Eq.  37,  this  can  be  written 

y  =  U''au'^x 
Or 

U^y  =  AU'^x  (38) 

Let  us  define  "canonical"  incident  and  reflected  waves,  y'  and  x'  by 


x'  =  U^x  ;  X  =  Ux' 

y'  =  U^y  ;  y  =  U%' 

(39) 

Then  also 

y'  =  Ax' 

(40a) 

and  in  components 

4  =  Vk 

(40b) 

But  what,  in  fact,  are  these  "canonical  waves"?  Well,  first  the  unitary  matrix  U 
in  Eq.  39  may  be  regarded  as  a  transfer  matrix  for  a  lossless  "black  box",  with 
M  ports  on  each  side,  such  that  putting  in  the«|^xj^j  at  one  end  gives  the^Xj^^  at 
the  other  (with  no  reflections),  and  vice-versa,  if  the  reflections  from  the 
antenna  y^^J  are  traveling  back  towards  this  box,  the  ^^yj^  emerges  without 
losses.  Obviously  such  a  "black  box"  satisfies  all  our  conservation  of  energy 
limitations,  and  we  should  not  be  surprised  to  find  this  degree  of  flexibility 
in  our  final  equations.  (One  can  dispute  whether  such  a  lossless,  reflectionless 
black  box  is  realizable;  our  equations  only  say  that  if  realizable,  it  can  be 
inserted. ) 

Further,  each  of  the  x^  represents  a  particular  weighted  sum  of  the  actual 
incident  waves 


x: 


-  2  Uu .X. 
;  kj  J 


(41) 
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and  inversely,  given  the  set^x^  , 

^  ^  (42) 

The  set  of  excitations,  >  corresponding  to  one  of  the  x'  non- vanishing  and  all 

others  zero,  obviously  defines  some  combination  of  beams,  i,  e,  ,  some  radiation 

pattern.  Thus,  each  of  the  x'  corresponds  to  exciting  some  particular  radiation 

(1)  .  ... 

pattern,  and  (referring  to  Part  I'  ')  it  is  readily  observed  that  since  the  weightings  in 

Eq,  42  are  from  a  unitary  matrix  which  diagonalizes  F,  the  set  of  corresponding 

"canonical”  radiation  patterns  (beam  combinations)  are  in  fact  orthogonal  to  each  other. 

Since  they  are  orthogonal,  their  "feed-lines"  are  decoupled;  as  shown  in  Eq,  40b,  a 

particular  xj  gives  rise  only  to  the  corresponding  "reflection",  yj  in  its  own  feed-line. 

The  further  arbitrariness  in  U  corresponds  to  the  statement  that  certain  "canonical" 

pattern  excitations  are  equivalent  as  far  as  the  reflections  produced  (equal  Xj^)t  and 

hence  as  far  as  efficiency  in  radiating  the  energy  incident  in  the  feeds.  Thus  the 

"canonical"  pattern  representation  is  not  unique;  in  particular,  even  if  all  the  were 

unequal,  phase  shifters  in  the  lines  (or  equivalently,  shift  in  choice  of  phase  references) 

would  in  no  manner  change  the  basic  energy  flow  relationships.  It  is  also  interesting  to 

note  that  the  eigenvalues  y  (of  F)  are  directly  interpretable  in  these  canonical  excitations, 

in  terms  of  the  reflection  coefficients,  through  the  relations 

\  =  y  >  -  I'k 

It  is  not  clear  to  what  extent  the  "canonical"  patterns  may  be  of  direct  physical 
interest.  Presumably,  our  physical  interest  lies  in  the  set  of  beam  patterns  which  we 
specified  at  the  outset.  Although  Eq.  42  also  can  be  regarded  as  describing  how  to 
achieve  this  set  by  simultaneous  excitation  at  a  number  of  terminals,  each  of  which  is 
decoupled  from  the  other,  the  totality  of  waves  for  a  single  desired  beam  involves  all 
terminals,  and  no  obvious  engineering  advantage  ensues.  What  appears  to  be  more  to 
the  point  is  a  study,  within  some  specific  systems  applications  of  multiple -beam 
antennas,  of  the  possibilities  for  finding  a  "minimum"  S  in  some  sense.  The  important 
result  is  that  the  character  of  S  is  largely  fixed  by  the  beam  design,  and  it  could  be 
important  to  an  antenna  design  engineer  to  understand  what  aspects  of  the  cross -couplings 
(elements  of  S)  cannot  be  changed  without  altering  the  beam  structure. 

6.  Need  for  Further  Study 

The  results  of  this  paper  can  only  be  regarded  as  a  first  step  in  character¬ 
izing  the  cross -couplings  within  multiple-beam  antenna  systems.  It  appears  that 
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more  explicit  relationships  ought  to  be  determined,  probably  in  relation  to 
specific  system  applications,  in  order  to  arrive  at  the  really  desirable  point  for 
the  antenna  design  engineer;  to  determine  what  limitations  he  faces  a  priori  in 
trying  to  modify  an  antenna  structure  so  as  to  reduce  cross-couplings  without 
exceeding  radiation  pattern  tolerances  specified  by  his  customer.  For  example, 
in  Part  I  it  was  pointed  out  that  much  of  the  cross -coupling  appears  to  be  controlled 
by  side-lobe  effects,  where  perhaps  more  freedom  is  possible.  The  important 
factor  would  be  the  detailed  nature  of  the  connection,  and  this  we  do  not  claim  to 
have  truly  elucidated  in  this  initial  formulation. 

What  this  paper  has  achieved,  we  believe,  is  to  make  explicit  the  basic 
mathematical  structure  which  can  be  a  starting  point  for  such  further  research. 
Over  and  above  such  additional  research  in  even  the  "simple"  lossless,  reciprocal 
case,  lie  further  explorations  pointed  out  earlier  in  connection  with  such  system 
properties  as  "lossy",  "non-reciprocal",  and  the  like. 
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APPENDIX  1:  Lack  of  Uniqueness  in  Unitary  Diagonalizing  Matrices 

Consider  an  MxM  Hermitian  matrix  F.  Its  eigenvectors  form  a  mutually 

orthogonal  set,  or,  for  a  degenerate  eigenvalue  can  be  constructed  as  such  (the 

Gram-Schmidt  procedure).  Further,  the  eigenvectors  can  be  normalized  to  unit 

length.  Thus,  if  denotes  the  orthonormalized  eigenvector  (with  components 

(k) 

,  m  =  1,  ....  M),  we  can  write  the  matrix  statement 
m 


^  (Al-1) 

(k), 

A  unitary  matrix  U  with  the  set  '•  as  its  columns  then  has  the  property  of 

diagonalizing  F  (Eq.  4) 

U'^FU  y  (Al-2) 

However,  any  eigenvector  can  still  be  multiplied  in  all  its  components  by  any 
complex  number  of  unit  magnitude,  i.  e.  ,  by  a  factor  of  the  form  exp  (iO),  without 
affecting  the  statement  in  Eq.  Al-1,  The  corresponding  effect  on  the  construction 
of  U  is  that  U  in  Eq.  Al-2  can  be  replaced  by  any  unitary  matrix,  V,  of  the  form 


V  =  UD 

where  D  is  a  diagonal  unitary  matrix  with  the  form 


(Al-3) 


/exp(i0^) 


D  = 


exp(i02) 


(Al-4) 


and  where  the  0^,  d^.  .  .  .  are  quite  a.rbitrary.  The  matrix  D  clearly  has  the 
properties  of  being  both  unitary  and  symmetric, 

T 

D  =  D 


D'^D  =  I 


(Al-5) 


where  I  is  the  identity  matrix,  whose  elements  are 

Next,  suppose  F  has  a  degenerate  (multiple)  eigenvalue.  If  the  order  of 
degeneracy  (multiplicity)  of  a  particular  eigenvalue  is  r,  then  the  secular  equation 


Al-1 
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determining  the  corresponding  eigenvectors  gives  additional  r-1  undetermined 
coefficients  (i.  e.  ,  in  addition  to  the  usual  indeterminancy  within  a  complex  scale 
factor  for  the  vector  as  a  whole).  One  is  then  free  to  choose  these  coefficients  to 
form  r  linearly  independent  eigenvectors  (automatically  orthogonal  to  the  eigenvectors 
corresponding  to  any  of  the  other  eigenvalues),  and  in  particular  usually  to  construct 
them  as  a  mutually  orthogonal  and  normalized  set  (e,  g.  ,  the  Gram-Schmidt 
procedure)  of  r  eigenvectors,  which  can  then  in  turn  be  used  in  constructing  U. 
However,  given  one  such  set,  any  other  set  of  r  linear  combinations  of  these, 
constructed  so  as  to  provide  orthonormality  will  also  be  suitable  for  constructing  U. 
An  orthonormality-preserving  traTisformation  from  one  subset  to  another  is,  of 
course,  just  an  arbitrary  unitary  transformation  in  the  subspace  spanned  by  the 
eigenvectors  of  the  particular  eigenvalue.  I.  e.  ,  in  the  case  of  degenerate 
eigenvalues,  we  must  extend  Eq.  Al-3  by  the  statement  that  more  generally,  given 
any  U,  we  can  equally  well  satisfy  Eq.  Al-2  by  any 


W  =  UD 


(Al-6) 


where  D  has  the  quasi -diagonal  partitioned  form 


x(l) 


.(2) 


D 


(3)1 


D  =1 


(Al-7) 


-[5(n^)-, 


The  partition, 


is  any  unitary  matrix  of  order  r.  X  r.,  where  r.  is  the 

.  th  J  J  J 

order  of  degeneracy  (multiplicity)  of  the  j  different  eigenvalue  of  F,  in  some 

assumed  ordering  of  the  eigenvalues.  In  this  form,  the  columns  of  include  the 

arbitrary  phase  factors  outlined  in  Eq.  Al-4;  in  fact  for  any  non-degenerate 

eigenvalue,  the  corresponding  partition  is  only  a  single  element  of  general 

form  exp 


It  is  to  be  noted  that  if  the  matrix  y  is  represented  in  a  partitioned  form 
similar  to  Eq.  Al-7,  the  partitions  are  a  set  of  identity  matrices,  multiplied 
by  scalars. 


Al-2 
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and  hence  y  commutes  with  any  D  matrix  of  the  form  given  in  Eq.  Al-7, 


yD  =  Dy 


(Al-8) 


(Al-9) 
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APPENDIX  2  General  Solution  for  S  in  Singular  Case 

If  some  of  the  in  Eq,  7  are  zero,  then  A  ^  is  strictly  undefined,  and  we 
cannot  proceed  as  in  Eq.  11a.  However,  we  can  certainly  still  write  Eq.  11b, 

S  =  TAU" 

to  introduce  a  new  matrix  T.  If  some  of  the  diagonal  elements  of  the  diagonal 
matrix  of  A  vanish,  then  it  is  readily  observed  that  in  the  matrix  multiplication  of 
Eq.  A2-1,  the  corresponding  columns  of  T  (whatever  their  values  maybe)  never 
contribute  to  the  elements  of  S;  likewise  the  corresponding  rows  of  (or  columns  of 
U)  never  enter.  This  is  really  the  implication  of  the  singularity  in  A,  that  the  new 
unknown  matrix  T  can  never  be  completely  specified  from  a  knowledge  of  S.  However, 
since  S  rather  than  T  is  our  true  objective,  this  is  no  hindrance  to  our  solution. 

Further,  as  in  Eq.  13,  we  can  also  introduce  a  new  vector  f, 

S  =  Az 

=  z^A'  (A2-2) 

in  which,  by  the  property  of  some  of  the  vanishing,  the  corresponding  elements 
of  5  always  vanish.  Thus  this  f  is  arbitrary  only  in  those  elements  corresponding 
to  non -vanishing  X^^. 

We  thus  again  derive  the  inequality  of  Eq.  14. 

S'S>f"T-TS  (A2-3) 

But  now,  due  to  the  zero  elements  just  mentioned  in  the  class  of  otherwise  arbitrary 
J-vectors,  this  inequality  (as  is  readily  observed  by  inspection)  in  fact  specifies 
limitations  on  all  columns  of  T,  except  exactly  those  already  mentioned  as  being 

4-  4-  x  x 

irrelevant  to  the  form  of  S.  The  matrix  T  T  is  always  Hermitian,  since  (T  T)  =  T  T, 
and  positive  definite  since  Z  T'TZ=(TZ)  (TZ)>0,  so  a  representation  as  in  Eq.  1  5 
is  always  possible 

t't  -  pV^P  (A2-4) 

2 

where  P  is  again  unitary,  and  tx  diagonal  with  elements  positive  real  or  zero.  But 
now  (inserting  a  step  previously  left  to  the  reader  between  Eqs.  15a  and  16), 

Eq.  A2-3  becomes 

S^f>£^PV^P?  (A2-5) 

We  can  now  introduce  the  transformation 


t?  =  Pf 


(A2-6) 
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where  77  can  be  any  arbitrary  vector 


P% 


lovitlcj.  that  in  the  inverse  transformation 


(A2-7) 


the  required  elements  of  f  always  vanish.  This  implies  that  for  each  arbitrary  77, 
the  appropriate  rows  of  (columns  of  P)  must  always  be  such  as  to  produce  this 
vanishing  of  the  required  elements  of  S.  But,  regarding  this  as  a  specification 
to  be  satisfied  by  any  P  providing  a  solution  through  Eq.  A2-4,  the  inequality  of 
Eq.  A2-5  becomes  just 


+  +2 

77  77  >  77  p  77 


(A2-8) 


Or,  again, 

^  1  for  all  k  (A2-9) 

We  can  certainly  again  factor  Eq.  A2-4,  to  give 


T  =  WpP 


where  W  is  any  unitary  matrix;  and  thus  from  Eq.  A2-1,  write  (exactly  as  in  Eq.  19) 

S  =  WmEAU"^  (A2-10) 

Finally,  we  observe  that  the  peculiar  limitation  on  certain  columns  of  P 
(especially  peoaliar -looking  now  since  it  is  related  to  the  form  of  an  arbitrary 
vector,  77,  which  does  not  enter  our  final  equations)  is  in  fact  irrelevant  since 
P  enters  only  through  the  product  PA,  and  it  is  exactly  these  columns  of  P  which 
do  not  contribute  to  the  computation  of  S  because  of  the  vanishing  of  the  corresponding 

V/e  have  thus  demonstrated  our  assertion  that  the  result,  Eq.  19,  holds 
without  change  in  the  limiting  case  where  A  is  singular  through  the  vanishing 
of  some  of  its  diagonal  elements. 
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2  copies 

M  6 

AFCRL  AFRD  (CRRELT) 

L,  G.  Hanscom  Field,  Bedford,  Mass. 

10  copies 

M  17 

AFCRL  AFRD  (CRRD) 

Attn:  Carlyle  J.  Sletten 

L.  G.  Hanscom  Field,  Bedford,  Mass. 

3  copies 

AFCRL  AFRD  (CRRD) 

Attn:  Contract  Files 

L.  G.  Hanscom  Field,  Bedford,  Mass. 

2  copies 

M  54 

Hq.  AFCCDD  (CCRRC,  Capt.  John  J.  Hobson) 

L.  G.  Hanscom  Field,  Bedford,  Mass. 

M  59 

AFCCDD 

L.  G,  Hanscom  Field,  Bedford,  Mass. 

N  1 

Director,  Avionics  Division  (AV) 

Bureau  of  Aeronautics 

Department  of  the  Navy 

Washington  25,  D.  C. 

N  3 

Chief,  Bureau  of  Shi  os 

Department  of  the  Navy 

Washington  25,  D.  C. 

Attn:  Code  690 

(3  ode 


N  9 


N  16 


N  23 


N  26 


N  27 


N  29 


N  30 


N  35 


N  37 


N  48 


N  85 


Organization 


Chief,  Bureau  of  Ordnance 
Department  of  the  Navy 
Washington  25,  D.  C. 

Attn:  Code  Ad3 

Commander 

U.  S.  Naval  Air  Missile  Test  Center 
Point  Mugu,  California 
Attn:  Code  366 

U.  S.  Naval  Ordnance  Laboratory 
White  Oak,  Silver  Spring  19,  Maryland 
Attn:  The  Library 

Commander 

U.  S.  Naval  Ordnance  Test  Station 
Cnina  Lake,  California 
Attn:  Code  753 

Librarian 

U.  S.  Naval  Postgraduate  School 
Monterey,  California 

Director 

U.S.  Naval  Research  Laboratory 
Washington  25,  D,  C, 

Attn:  Code  2027 

Dr,  J,  I,  Bohnert,  Code  5210 
U.S.  Naval  Research  Laboratory 
Washington  25,  D,  C, 

Commanding  Officer  and  Director 

U.  S.  Navy  Underwater  Sound  Laboratory 

Fort  Trumbull,  New  London,  Connecticut 

Chief  of  Naval  Research 
Department  of  the  Navy 
Washington  25,  D.  C. 

Attn;  Code  427 

Commanding  Officer 
U.  S.  Naval  Air  Development  Center 
Johnsville,  Pennsylvania 
Attn:  NADC  Library 

Commanding  Officer  and  Director 

U.S.  Navy  Electronics  Laboratory  (Library) 

San  Diego  52,  California 


2  copies 


5 


Code 
N  86 


N  87 

N  88 

N  91 

N  92 

N  93 

N  123 

N  141 

N  73 

I  601 


Organization 


Chief,  Bureau  of  Ordnance 
Department  of  the  Navy 
Surface  Guided  Missile  Branch 
Washington  25,  D.  C. 

Attn:  Code  ReSl 

Chief,  Bureau  of  Ordnance 
Department  of  the  Navy 
Washington  25,  D.  C. 

Attn;  Fire  Control  Branch  (ReS4) 

Department  of  the  Navy 
Bureau  of  Aeronautics 
Technical  Data  Division,  Code  4106 
Washington  25,  D.  C. 

Commander 

U.  S.  Naval  Air  Test  Center 
Patuxent  River,  Maryland 
Attn:  ET-315,  Antenna  Branch 

New  York  Naval  Shipyard 
Material  Laboratory,  Code  932 
Brooklyn  1,  New  York 
Attn;  Mr.  Douglas  First 

Commanding  Officer 

U.  S.  Naval  Ordnance  Laboratory 

Corona,  California 

Attn;  Documents  Librarian 

Chief,  Bureau  of  Ships 
Department  of  the  Navy 
Washington  25,  D.  C. 

Attn:  Code  817B 

ARDC  Regional  Office 
c/o  Department  of  the  Navy 
Room  4549,  Munitions  Building 
Washington  25,  D.  C. 

Office  of  Naval  Research  10  copies 

Branch  Office,  London 

Navy  100,  Box  39 

F.P.  O.  New  York,  N.  Y. 

Aero  Geo  Astro  Corp. 

1200  Duke  Street 
Alexandria,  Virginia 
Attn:  Library 


6 


Code 
I  1 

I  388 

I  3 

I  374 

I  205 

I  8 

I  469 

I  13 

I  246 

I  940 


Organization 


Airborne  Instruments  Laboratory,  Inc. 

Division  of  Cutler  Hammer 
Walt  Whitman  Road 
Melville,  L.  I,  New  York 
Attn:  Library 

Air  com,  Inc. 

354  Main  Street 
Winthrop,  Mass. 

Andrew  Alfor<i  Consulting  Engineers 
299  Atlantic  Avenue 
Boston  10,  Mass. 

ACF  Industries,  Inc. 

Southeast  Corner,  52nd  Avenue  and  Jackson  Street 
Bladensburg,  Prince  George's  County,  Maryland 
Attn:  Librarian 

Battelle  Memorial  Institute 
505  King  Avenue 
Columbus  1,  Ohio 

Attn:  Wayne  E.  Rife,  Project  Leader 
Electrical  Engineering  Division 

Bell  Aircraft  Corporation 

Post  Office  Box  One 

Buffalo  5,  New  York 

Attn:  Eunice  P,  Hazelton,  Librarian 

Bell  Telephone  Laboratories 
Murray  Hill 
New  Jersey 

Bell  Telephone  Laboratories,  Inc, 

Whippany  Laboratory 

Whippany,  New  Jersey 

Attn:  Technical  Information  Library 

Bendix  Corporation 
Pacific  Division 
11600  Sherman  Way 
North  Hollywood,  California 
Attn:  Engineering  Library 

Aerospace  Corp. 

Box  95085 

Los  Angeles  45,  California 
Attn:  Library 


Code 
I  247 


I  248 

I  249 

I  250 

I  892 

I  252 

I  253 

I  470 

I  918 

I  919 


Organization 


Bendix  Radio  Division 

Bendix  Aviation  Corporation 

E.  Joppa  Road 

Tow  son  4,  Maryland 

Attn:  Dr.  D.  M.  Allison,  Jr. 

Director  Engineering  &  Research 

Bjorksten  Research  Laboratories,  Inc. 

P.  O.  Box  265 
Madison,  Wisconsin 
Attn:  Librarian 

Boeing  Airplane  Company 
Pilotless  Aircraft  Division 
P.O.  Box  3707 
Seattle  24,  Washington 

Attn:  R.  R.  Barber,  Library  Supervisor 

Boeing  Airplane  Company 

Wichita  Division  Engineering  Library 

Wichita  1,  Kansas 

Attn:  Kenneth  C.  Knight,  Library  Supervisor 

Boeing  .\irplane  Company 
7755  East  Marginal  Way 
Seattle  24,  Washington 
Attn:  Mr,  D.  L.  Scidmore 

Brush  Beryllium  Company 
4301  Perkins  Avenue 
Cleveland  3,  Ohio 
Attn:  N.  W.  Bass 

Chance  Vought  Aircraft,  Inc. 

9314  West  Jefferson  Street 
Dallas,  Texas 

Attn:  A.  D.  Pattullo,  Librarian 

Chu  Associates 
P.  O.  Box  387 
Whitcomb  Avenue 
Littleton,  Mass. 

Collins  Radio  Co. 

855  35th  Street,  N.  E. 

Cedar  Rapids,  Iowa 
Attn:  Dr.  R.L.  McCreary 

Collins  Radio  Co. 

1200  North  Alma  Road 
Richardson,  Texas 
Attn:  C.  D.  Tipton 


2  copies 


Code 
I  126 

I  254 

I  25 

I  255 

I  28 

I  257 

I  258 

I  259 

I  187 

I  415 


Organization 


Convair,  A  Division  of  General  Dynamics  Corp. 
Fort  Worth,  Texas 
Attn;  K.  G.  Brown 

Division  Research  Librarian 

Convair,  A  Division  of  General  Dynamics  Corp. 
3165  Pacific  Highway 
San  Diego  12,  California 
Attn;  Ivirs.  Dora  B.  Burke, 

Engineering  Librarian 

Cornell  Aeronautical  Laboratory,  Inc. 

4455  Genesee  Street 
Buffalo  21,  New  York 
Attn;  Librarian 

Dalmo  Victor  Company 
A  Division  of  Textron,  Inc. 

1515  Industrial  Way 
Belmont,  California 

Attn;  Mary  Ellen  Addems,  Technical  Librarian 

Dome  and  Margolin,  Inc. 

29  New  York  Avenue 
Westbury,  Long  Island,  N.  Y. 

Douglas  Aircraft  Co.  ,  Inc. 

827  Lapham  Street 
El  Segundo,  California 
Attn;  Engineering  Library 

Douglas  Aircraft  Company,  Inc, 

3000  Ocean  Park  Boulevard 
Santa  Monica,  California 
Attn:  Peter  Duyan,  Jr. 

Chief,  Electrical/ Electronics  Section 

Douglas  Aircraft  Company,  Inc. 

2000  North  Memorial  Drive 
Tulsa,  Oklahoma 

Attn;  Engineering  Librarian,  D-250 

Electromagnetic  Research  Corporation 
5001  College  Avenue 
College  Park,  Maryland 
Attn;  Mr.  Martin  Katzin 

Electronics  Communication 
1830  York  Road 
Timonium,  Maryland 


Code 
I  299 

I  204 

I  262 

I  147 


I  264 

I  266 

I  269 

I  793 

1  893 


Organization 


Electronic  Specialty  Company 
5121  San  Fernando  Road 
Los  Angeles  39,  California 
Attn;  Donald  L.  Ma.rgerum 

Chief  Engineer,  Radiating  Systems  Division 

Emerson  and  Cuming,  Inc, 

869  Washington  Street 
Canton,  Mass, 

Attn;  Mr.  W,  Cuming 

Emerson  Electric  Mfg,  Co, 

8100  West  Florissant  Avenue 

St.  Louis  21,  Missouri 

Attn;  Mr,  E.R,  Breslin,  Librarian 

Emerson  Radio-Phonograph  Corp, 

Emerson  Research  Laboratories 
701  Lamont  Street,  N.  W. 

Washington  10,  D.  C. 

Attn:  Mrs.  R,  Corbin,  Librarian 

Classified  to:  Mrs.  Coleen  G,  Molineau 

Fairchild  Aircraft-Missiles  Division 
Fairchild  Eng,  and  Airplane  Corp. 

Hagerstown,  Maryland 
Attn:  Library 

Federal  Telecommunication  Laboratories 
500  Washington  Avenue 
Nutley  10,  New  Jersey 
Attn;  Technical  Library 

General  Electric  Company 
Electronics  Park 
Syracuse,  New  York 
Attn:  Documents  Library 

B.  Fletcher,  Building  3-143A 

General  Electric  Company 

Missile  and  Space  Vehicle  Department 

3198  Chestnut  Street,  Philadelphia,  Pennsylvania 

Attn:  Documents  Library 

General  Electric  Company 
3750  D  Street 

Philadelphia  24,  Pennsylvania 
Attn:  Mr.  H.G.  Low 

Missile  and  Space  Vehicle  Department 


Code 
I  270 

I  48 

I  448 

I  272 

I  273 

I  737 

I  311 

I  274 

I  207 

I  56 


Organization 


General  Precision  Laboratory,  Inc. 

63  Bedford  Road 
Pleasantville,  New  York 
Attn:  Librarian 

Goodyear  Aircrjift  Corp. 

1210  Massillon  Road 

Akron  15,  Ohio 

Attn:  Library,  Plant  G 

Granger  Associates 

Electronic  Systems 

974  Commercial  Street 

Palo  Alto,  California 

Attn:  John  V.  N.  Granger,  President 

Grumman  Aircrzift  Engineering  Corporation 

Bethpage,  Long  Island,  New  York 

Attn:  Engineering  Librarian,  Plant  No,  5 

Hallicrafter s  Company 

4401  West  5th  Avenue 

Chicago  24,  Illinois 

Attn:  La  Verne  LaGioia,  Librarian 

The  Hallicrafters  Co. 

5th  and  Kostner  Avenues 
Chicago  24,  Illinois 

Attn:  Henri  Hodara,  Head  of  Space  Connmunication 

Hermes  Electronics  Company 
75  Cambridge  Parkway 
Cambridge,  Mass. 

Attn:  Mrs.  Lois  Seulowitz 
Technical  Librarian 

Hoffman  Electronics  Corp, 

37  61  South  Hill  Street 
Los  Angeles  7,  California 
Attn:  Engineering  Library 

Hughes  Aircraft  Company 
Antenna  Department 
Building  12,  Maiil  Station  2714 
Culver  City,  California 
Attn:  Dr.  W.  H.  Kummer 

Hughes  Aircraft  Company 
Florence  and  Teale  Streets 
Culver  City,  California 
Attn:  Mr.  Louis  L.  Bailin 

Manager,  Antenna  Department 
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Code 
I  302 

I  414 

I  265 

I  230 

I  241 

I  279 


I  468 

I  136 

I  280 

I  63 


Organization 

International  Business  Machines  Corp. 

Space  Guidance  Center- Federal  Systems  Division 
Owego,  Tioga  County,  New  York 
Attn:  Technical  Reports  Center 

International  Resistance  Company 
401  N.  Broad  Street 
Philadelphia  8,  Pa. 

Attn:  Research  Library 

ITT  Laboratories 
3700  East  Pontiac  Street 
Fort  Wayne  1,  Indiana 
Attn:  Technical  Library 

Jansky  and  Bailey,  Inc. 

1339  Wisconsin  Avenue,  N.  W. 

Washington  7,  D.  C. 

Attn:  Mr.  Delmer  C.  Ports 

Dr.  Henry  Jasik,  Consulting  Engineer 
298  Shames  Drive 
Brush  Hollow  Industrial  Park 
Westbury,  New  York 

Lockheed  Aircraft  Corporation 

2555  N.  Hollywood  Way 

California  Division  Engineering  Library 

Department  72-25,  Plant  A-1,  Building  63-1 

Burbank,  California 

Attn:  N.  C.  Harnois 

Lockheed  Aircraft  Corporation 

Missile  Systems  Division  Research  Library 

Box  504,  Sunnyvale,  California 

Attn;  Miss  Eva  Lou  Robertson,  Chief  Librarian 

The  Martin  Company 

P.  O,  Box  179 

Denver  1,  Colorado 

Attn;  Mr.  Jack  McCormick 

The  Martin  Company 
Baltimore  3,  Maryland 
Attn;  Engineering  Library 
Antenna  Design  Group 

Mathematical  Reviews 
190  Hope  Street 
Providence  6,  Rhode  Island 


Code 
I  66 

I  282 

I  283 

I  116 

I  471 

I  390 

I  379 

I  145 

I  648 

I  934 


Organization 

The  W.L.  Maxson  Corporation 
47  5  10th  Avenue 
New  York,  New  York 
Attn:  Miss  Dorothy  Clark 

McDonnell  Aircraft  Corporation,  Dept,  644 

Box  516,  St.  Louis  66,  klissouri 

Attn:  C.  E.  Zoller,  Engineering  Library 

McMillan  Laboratory,  Inc, 

Brownville  Avenue 
Ipswich,  Mass. 

Attn:  Security  Officer,  Document  Room 
Melpar,  Inc. 

3000  Arlington  Boulevard 

Falls  Church,  Virginia 

Attn:  Engineering  Technical  Library 

Microwave  Associates,  Inc. 

South  Avenue 
Burlington,  Mass. 

klicrowave  Development  Laboratory 
90  Broad  Street 
Wellesley  57,  Mass. 

Attn:  N.  Tucker,  General  Manager 

Microwave  Engineering  Laboratories,  Inc. 

943  Industrial  Avenue 

Palo  Alto,  California 

Attn:  J.  L.  Melchor,  President 

Microwave  Radiation  Company,  Inc. 

11121  Hindry  Avenue 

Los  Angeles  45,  California 

Attn:  Mr.  Morris  J,  Ehrlich,  President 

'T’he  Mitre  Corporation 
244  Wood  Street 
Lexington  73,  Mass, 

Attn:  Mrs.  Jean  E.  Claflin,  Librarian 
Motorola,  Inr.. 

Phoenix  Research  Laboratory 
3102  N,  56th  Street 
Phoenix,  Arizona 
Attn:  Dr.  A.  L.  Aden 


Code 
I  641 

I  284 


I  285 

I  82 

I  286 

1  287 


I  225 

I  288 

I  232 

I  289 


Organization 


National  Research  Council 

Radio  &  Electrical  Engineering  Division 

Ottawa,  Ontario,  Canada 

Attn;  Dr.  G.  A,  Miller,  Head  Microwave  Section 

North  American  Aviation,  Inc. 

12214  Lakewood  Boulevard 
Downey,  California 

Attn;  Technical  Information  Center  (495-12) 
Space  &  Information  Systems  Division 

North  American  Avaiation,  Inc, 

Los  Angeles  International  Airport 
Los  Angeles  45,  California 
Attn;  Engineering  Technical  File 

Northrop  Corporation 
Nor  air  Division 
1001  E.  Broadway 
Hawthorne,  California 

Attn;  Mr,  E.  A,  Freitas,  Library  Dept,  3145 

Page  Communications  Engineers,  Inc, 

2001  Wisconsin  Avenue,  N.  W. 

Washington  7,  D.  C. 

Attn:  Mrs.  Ruth  Temple,  Librarian 

Philco  Corporation 
C  and  Tioga  Streets 
Philadelphia  34,  Pennsylvania 
Attn:  Mrs.  Dorothy  B.  Collins 
Research  Librarian 

Pickard  and  Burns,  Inc. 

240  Highland  Avenue 
Needham  94,  Mass. 

Attn:  Dr.  Richard  H.  Woodward 

Polytechnic  Research  &  Development  Co,  Inc. 

202  Tillary  Street 
Brooklyn  1,  New  York 
Attn;  Technical  Library 

Radiation  Engineering  Laboratory 
Main  Street 
Maynard,  Mass. 

Attn:  Dr,  John  Ruze 

Radiation,  Inc. 

P.  O,  Drawer  37 
Melbourne,  Florida 

Attn:  Technical  Library,  Mr.  M.  L,  Cox 


Code 
I  914 

I  233 

I  290 

I  291 


I  473 

I  757 

I  789 

I  930 

I  292 


Organization 

Radiation  Systems,  Inc. 

440  Swann  Avenue 
Alexandria,  Virginia 
Attn:  Library 

Radio  Corporation  of  America 
RCA  Laboratories 
Rocky  Point,  New  York 
Attn:  P.  S.  Carter,  Lab.  Library 

RCA  Laboratories 

David  Sarnoff  Research  Center 

Princeton,  New  Jersey 

Attn:  Miss  Fern  Cloak,  Librarian 

Radio  Corporation  of  America 
Defense  Electronic  Products 
Building  10,  Floor  7 
Camden  2,  New  Jersey 

Attn:  Mr.  Harold  J.  Schrader,  Staff  Engineer 

Organization  of  Chief  Technical  Administrator 

Radio  Corporation  of  America 
Missile  Control  &  Electronics  Division 
Bedford  Street,  Burlington,  Mass. 

Attn:  Librarian 

Radio  Corporation  of  America 
Surface  Communications  Systems  Laboratory 
75  Varick  Street 
New  York  13,  N.  Y. 

Attn;  Mr.  S.  Krevsky 

Radio  Corporation  of  America 

West  Coast  Missile  and  Surface  Radar  Division 

Engineering  Library,  Building  306/2 

Attn:  L.  R.  Hund,  Librarian 

8500  Balboa  Boulevard,  Van  Nuys,  Ccdifornia 

Radio  Corp.  of  .iXmerica 
Defense  Electronic  Products 
Advanced  Military  Systems 
Princeton,  New  Jersey 
Attn;  Mr,  David  Shore 

Director,  USAF  Project  RAND 

Via;  AF  Liaison  Office 

The  Rand  Corporation 

1700  Main  Street 

Santa  Monica,  California 


Code 
I  547 

I  373 

I  293 

I  294 

I  472 

I  510 

I  295 

I  184 

I  391 

I  142 


Organization 

The  Rand  Corporation 
1700  Main  Street 
Santa  Monica,  California 
Attn:  Technical  Library 

Rantec  Corporation 

23999  Ventura  Boulevard 

Calabasas,  California 

Attn:  Grace  Keener,  Office  Manager 

Raytheon  Company 

State  Road,  Wayland  Laboratory 

Wayland,  Mass. 

Attn:  Mr.  Robert  Borts 

Raytheon  Company 
Wayland  Laboratory 
Wayland,  Mass. 

Attn:  Miss  Alice  G.  Anderson,  Librarian 

Raytheon  Company 
Missile  Systems  Division 
Hartwell  Road,  Bedford,  Mass, 

Attn:  Donald  H.  Archer 

Remington  Rand  Univ. 

Division  of  Sperry  Rand  Corporation 
1900  West  Allegheny  Avenue 
Philadelphia  29,  Pa. 

Attn:  Kto.  John  F.  McCarthy 

R  &  D  Sales  and  Contracts 

Republic  Aviation  Corporation 
Farmingdale,  Long  Island,  New  York 
Attn:  Engineering  Library 

Ryan  Aeronautical  Company 
2701  Harbor  Drive 
Lindbergh  Field 
San  Diego  12,  California 
Attn:  Library 

Sage  Laboratories,  Inc. 

3  Huron  Drive 
Natick,  Mass. 

Sanders  Associates,  Inc. 

95  Canal  Street 

Nashua,  New  Hampshire 

Attn:  Mr.  Norman  R.  Wild 


Code 
I  96 

I  546 

I  682 

I  312 

I  297 

I  367 

I  104 

I  260 

I  818 

I  240 

I  338 


Organization 

Sandia  Corporation,  Sandia  Base 

P.  O.  Box  5800,  Albuquerque,  New  Mexico 

Attn:  Classified  Document  Division 

Sandia  Corporation 

Attn:  Organization  1423 

Sandia  Base,  Albuquerque,  New  Mexico 

Scanwell  Laboratories,  Inc. 

6601  Scanwell  Lane 
Springfield,  Va. 

Space  Technology  Laboratories,  Inc. 

P.  O.  Box  95001 
Los  Angeles  45,  California 
Attn:  Mrs.  Margaret  C.  Whitnah 
Chief,  Information  Services 

Sperry  Gyroscope  Company 
Great  Neck,  Long  Island,  New  York 
Attn:  Florence  W.  Turnbull 
Engineering  Librarian 

Stanford  Research  Institute 
Documents  Center 
Menlo  Park,  California 
Attn:  Acquisitions 

Sylvania  Electric  Products,  Inc. 

100  First  Avenue,  Waltham  54,  Mass, 

Attn:  Charles  A.  Thornhill,  Report  Librarian 
Waltham  Laboratories  Library 

Sylvania  Electric  Products,  Inc. 

Electronic  Defense  Laboratory 
P.  O.  Box  205 
Mountain  View,  California 
Attn:  Library 

Sylvania  Reconnaissance  Systems  Laboratory 
Box  188,  Mountain  View,  California 
Attn:  Marvin  D.  Waldman 

TRG,  Inc. 

2  Aerial  Way 

Syosset,  New  York 

Attn:  M,  L.  Henderson,  Librarian 

A.  S.  Thomas,  Inc. 

355  Providence  Highway 
Westwood,  Mass. 

Attn:  A.  S.  Thomas,  President 
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Code 
I  708 

I  464 

I  139 

I  300 

U  1 

U  61 

U  157 

U  99 

U  3 


Organization 


Texas  Instruments,  Inc. 

6000  Lemmon  Avenue 
Dallas  9,  Texas 
Attn:  John  B.  Travis 

Systems  Planning  Branch 

Trans -Tech,  Inc. 

P.  O.  Box  346 

Rockville,  Maryland 

Attn:  Mr.  A.  C.  Blankenship 

Chief  Electronics  Engineer 

Westinghouse  Electric  Corp. 

Electronics  Division 
Friendship  Int'l  Airport  Box  1897 
Baltimore  3,  Maryland 
Attn:  Engineering  Library 

Zenith  Plastics  Company 
1600  W.  135th  Street 
Gardena,  California 
Attn:  Mr.  S.  S.  Oleesky 

Library  Geophysical  Institute  of  the 
Urdversity  of  Alaska 
College,  Alaska 

Brown  University 

Department  of  Electrical  Engineering 
Providence,  R.  I, 

Attn:  Dr,  C,  M,  Angulo 

California  Institute  of  Technology 
Jet  Propulsion  Laboratory 
4800  Oak  Grove  Drive 
Pasadena,  California 
Attn;  Mr.  I.  E.  Newlan 

California  Institute  of  Technology 
1201  E.  California  Street 
Pasadena,  California 
Attn:  Dr.  C.  Papas 

Space  Sciences  Laboratory 
Leuschner  Observatory 
University  of  California 
Berkeley  4,  California 

Attn:  Dr.  Samuel  Silver,  Professor  of  Engineering  Science 
and  Director,  Space  Sciences  Laboratory 
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Code 


U  100 


U  289 


U  239 


U  183 


U  238 


U  10 


U  86 


U  59 


U  102 


Organization 

University  of  California 
Electronics  Research  Laboratory 
332  Cory  Hall,  Berkeley  4,  California 
Attn:  J.R,  Whinnery 

University  of  Southern  California 
University  Park 
Los  Angeles,  California 
Attn:  Dr.  Raymond  L.  Chuan 

Director,  Engineering  Center 

Case  Institute  of  Technology 
Electrical  Engineering  Department 
10900  Euclid  Avenue 
Cleveland,  Ohio 

Attn:  Professor  Robert  Plonsey 

Columbia  University 

Department  of  Electrical  Engineering 

Morningside  Heights,  New  York,  N.  Y. 

Attn:  Dr.  Schlesinger 

University  of  Southern  California 
University  Park 
Los  Angeles  7,  California 
Attn:  7u.  K,  Kaprielian 

Associate  Professor  of  Electrical  Engineering 

Cornell  University 

School  of  Electrical  Engineering 

Ithaca,  New  York 

Attn:  Prof.  G.  C.  Dalman 

University  of  Florida 

Department  of  Electrical  Engineering 

Gainesville,  Florida 

Attn;  Prof.  M.  M.  Latour,  Library 

Library 

Georgia  Technology  Research  Institute 
Engineering  Experiment  Station 
722  Cherry  Street,  N.  W. 

Atlanta,  Georgia 

Attn;  Ml'S.  J.  H.  Crosland,  Librarian 

Harvard  University 
Technical  Reports  Collection 
Gordon  McKay  Library 
303  Pierce  Hall 

Oxford  Street,  Cambridge  38,  Mass. 

Attn:  Librarian 
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Code 
U  54 

U  103 

U  104 

U  169 

U  240 


U  22 

U  105 

U  228 

U  68 

U  32 


Organization 


Harvard  College  Observatory 
60  Garden  Street 
Cambridge  39.  Mass. 

Attn:  Dr.  Fred  L.  Whipple 

University  of  Illinois 
Documents  Division  Library 
Urbana,  Illinois 

University  of  Illinois 
College  of  Engineering 
Urbana,  Illinois 

Attn:  Dr,  P.  E,  Mayes,  Department  of  Electrical  Engineering 

Illinois  Institute  of  Technology 
3301  S.  Dearborn  Street 
Chicago  16,  Illinois 
Attn:  Dr.  George  I,  Cohn 

Illinois  Institute  of  Technology 
Technology  Center 

Department  of  Electrical  Engineering 
Chicago  16,  Illinois 
Attn:  Paul  C.  Yuen 

Electronics  Research  Laboratory 

The  Johns  Hopkins  University 
Homewood  Campus 
Baltimore  18,  l^ryland 

Attn:  Dr.  Donald  E.  Kerr,  Dept,  of  Physics 

The  Johns  Hopkins  University 
Applied  Physics  Laboratory 
8621  Georgia  Avenue 
Silver  Spring,  Maryland 
Attn:  Mr.  George  L.  Seielstad. 

University  of  Kansas 
Electrical  Engineering  Department 
Lawrence,  Kansas 
Attn:  Dr.  H,  Uns 

Lowell  Technological  Institute 

Research  Foundation 

P.  O.  Box  709,  Lowell,  Mass. 

Attn:  Dr.  Charles  R.  Mingins 

Massachusetts  Institute  of  Technology 
Research  Laboratory  of  Electronics 
Building  26,  Room  327 
Cambridge  39,  Mass. 

Attn:  John  H.  Hewitt 
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Code 
U  26 

U  34 

U  107 


U  79 

U  37 


U  108 

U  194 

U  39 

U  96 


Organization 


Massachusetts  Institute  of  Technology 
Lincoln  Laboratory 
P.  O.  Box  73 
Lexington  73,  Mass, 

Attn:  Mary  A,  Granese,  Librarian 

McGill  University 
Montreal,  Canada 
Attn:  Prof,  G,  A.  Woonton 

Director,  The  Eaton  Electronics 
Research  Laboratory 

University  of  Michigan 
Electronic  Defense  Group 
Engineering  Research  Institute 
Ann  Arbor,  Michigan 
Attn:  J,A,  Boyd,  Supervisor 

University  of  Michigan 
Engineering  Research  Institute 
Radiation  Laboratory 

912  N,  Main  Street,  Ann  Arbor,  Michigan 
Attn:  Prof,  K,  M,  Siegel 

University  of  Michigan 

Engineering  Research  Institute 

AVillow  Run  Laboratories,  Willow  Run  Airport 

Ypsilanti,  Michigan 

Attn:  Librarian 

University  of  Minnesota 

Minneapolis  14,  Minnesota 

Attn:  Mr,  Robert  H,  Stumm,  Library 

Physical  Science  Laboratory 

New  Mexico  College  of  Agriculture  and  Mechanic  Arts 
State  College,  New  Mexico 
Attn:  Mr,  H,W,  Haas 

New  York  University 

Institute  of  Mathematical  Sciences 

Room  802,  25  Waverly  Place 

New  York  3,  New  York 

Attn:  Professor  Morris  Kline 

Northwestern  University 
Microwave  Laboratories 
Evanston,  Illinois 
Attn:  R.  E.  Bean 
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U  78 

U  109 

U  45 

U  185 

U  48 

U  97 

U  184 

U  176 

u  no 

U  309 


Organization 

Ohio  State  University  Research  Foundation 
1314  Kinnear  Road 
Columbus  8,  Ohio 
Attn:  Dr.  T.  E.  Tice 

Department  of  Electrical  Engineering 

The  University  of  Oklahoma 
Research  Institute 
Norman,  Oklahoma 

Attn:  Prof.  C.  L.  Farrar,  Chairman 
Electrical  Engineering 

The  Pennsylvania  State  University 
Department  of  Electrical  Engineering 
University  Park,  Pennsylvania 

University  of  Pennsylvania 
Institute  of  Cooperative  Research 
3400  Walnut  Street,  Philadelphia,  Pa. 

Attn:  Department  of  Electrical  Engineering 

Polytechnic  Institute  of  Brooklyn 
Microwave  Research  Institute 
55  Johnson  Street,  Brooklyn,  New  York 
Attn:  Dr.  Arthur  A.  Oliner 

Polytechnic  Institue  of  Brooklyn 
Microwave  Research  Institute 
55  Johnson  Street 
Brooklyn,  New  York 
Attn:  Mr.  A.  E.  Laemmel 

Purdue  University 

Department  of  Electrical  Engineering 
Lafayette,  Indiana 
Attn;  Dr.  Schultz 

Stanford  University 

W.  W.  Hanser  Laboratory  of  Physics 
Stanford,  California 
Attn;  Microwave  Library 

Syracuse  University  Research  Institute 
Collendale  Campus 
Syracuse  10,  N.  Y. 

Attn:  Dr.  C.S.  Grove,  Jr., 

Director  of  Engineering  Research 

Technical  University 
Oestervoldgade  10  G 
Copenhagen,  Denmark 
Attn;  Prof.  Hans  Lottrup  Knudsen 
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Code 
U  186 

U  111 


U  51 

U  132 

U  133 

U  187 


Organization 


University  of  Tennessee 
Ferris  Hall 

W,  Cumberland  Avenue 
Knoxville  16,  Tennsesee 

The  University  of  Texas 

Electrical  Engineering  Research  Lab, 

P.  O.  Box  8026,  University  Station 
Austin  12,  Texas 
Attn:  Mr.  John  R,  Gerhardt 
Assistant  Director 

The  University  of  Texas 
Defense  Research  Laboratory 
Austin,  Texas 

Attn:  Claude  W.  Horton,  Physics  Librar 

University  of  Toronto 
Department  of  Electrical  Engr. 

Toronto,  Canada 
Attn:  Prof.  G.  Sinclair 

University  of  Washington 
Department  of  Electrical  Engineering 
Seattle  5,  Washington 
Attn:  G.  Held,  Associate  Professor 

University  of  Wisconsin 
Department  of  Electrical  Engineering 
Madison,  Wisconsin 
Attn:  Dr.  Scheibe 

U.  S.  Naval  Research  Lab. 

Washington  25,  D.  C. 

Attn:  A.  B.  Marston 

U.  S.  Army  Signal  Res.  and  Dev.  Labs. 
Belmar,  New  Jersey 
Attn;  A.  J.  Di Giacomo 

U.S.  Army  Signal  Res.  and  Dev.  Labs. 
Belmar,  New  Jersey 
Attn:  Boaz  Gelernter 

U.S.  Navy  Electronics  Lab. 

San  Diego  52,  California 
Attn:  B.  I.  Small 

RADC  (Attn;  A.  Feiner) 

Griffis s  AFB,  N.  Y. 


RADC  (Attn;  M.  Bickelhaupt) 
Griffis  s  AFB,  N.  Y. 


